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ABSTRACT 


In this survey and analysis of the meteoroid hazard to space 
vehicles in near-earth orbits," operations research methods have been 
used to establish and to propagate confidence levels with respect to all 
of the pertinent functional relations involving design and operations 
parameters . The graphically illustrated results show that the wall 
thickness, necessary for specified operational results, must be 
increased by a factor 2.00 to increase confidence from 50% to 75%, 
and that this factor would still be 1.81 if there were no uncertainty 
with respect to hypervelocity impact phenomena. 

Both the median values and the upper and lower fiducial limits for 
the 50% confidence interval for the mean flux, closing velocity relative 
to the earth, closing velocity relative to a vehicle in a near-earth 
orbit, density (specific gravity), just-puncturable thickness of hard 
aluminum, and just-puncturable thickness of hard stainless steel are 
shown for values of meteoroid mass between 10 grams and 
100 micro-micro grams. 


*For this report, near-earth orbits are considered to be of less 
than 500 km altitude. 
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Similar limits, the first, second, and third quartiles, are shown 
for puncture fluxes for walls of hard aluminum and hard stainless steel 
for thicknesses between 4 centimeters and 4 microns. 

Wall -thickness quartiles for 0.85, 0.90, and 0.95 probabilities of 
no puncture are shown for the hard aluminum and hard stainless steel 
walls with exposures between 1 square meter second and 1 mega-mega 
square meter second. s 

The reader should be aware that there still exists some 
controversy over the available meteoroid flux, density, velocity, 
and angular distribution statistics. Therefore, for any design 
studies the latest accepted values should be employed. 
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SUMMARY 


When G is a continuous function of the statistically independent 
chance variables x 1 , . . . , x n , then G is a random variable with the 
following sufficiently approximate standard deviation: 



The confidence C in the antilogarithm of an approximately normally 
distributed chance variable y is the same as the confidence in the 
variable; i. e . , 


y y - 0 . 6745cr y 

(C, 10 ) = (0.25, 10 Y ), (0.50, 10 ), 


(0. 75, 10 


y + 0. 6745 o'. 


y- 


)■ 


(29) 


*For this report, near-earth orbits are considered to be of less 
than 500 km altitude. 



For either sporadic or shower meteoroids, impact probability is 
sufficiently accurately approximated by either a Poisson or a piece- 
wise-Poisson description with respect to time, depending on whether 
the available information is specific or general, and on whether one 
is concerned with specific or average circumstances. But the 
combined effect for two such hazards also has a Poisson description: 

R = e“ A [^ u + ^ 12 + W*21 + <t >22 )t 2 + • • ■ + (4*ni + <>n 2 H n ] ^ 

The number of meteoroids with masses equal to or greater than 
m grams which should impact per second per square meter of hemi- 
spherical exposed area of a vehicle in a near-earth orbit is 


F > = 


10 


yi -i.i9 
1 m 


(37) 


where y x is an approximately normally distributed chance variable 
representing the uncertainty in the information and with mean yq and 
standard deviation <Xy as follows: 


y 1 = -12.86 


(44) 


Oy = 1.10 


(49) 


This mass dependence of flux is illustrated graphically in Fig. 1. 


The velocity in kilometers per second for a meteoroid of mass m 
grams, just before entering the earth' atmosphere, is 


v 


a 


10 y 2m°“* 


(56) 


where y 2 is an approximately normally distributed chance variable, 
which represents both a velocity spectrum and uncertainty in the 
information about it, with 


y 2 = 1.51 


( 61 ) 


o- = 0.12 
y 2 


(57) 


This mass dependence of velocity is illustrated graphically in Fig. 2. 
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Similarly, the closing velocity with respect to a vehicle in a 


earth orbit is 


c = 10 y 3m°' 041 

(62) 

3 = i' 52 

(84) 

< 

li 

o 
1 — 1 

(86) 


This mass dependence of closing velocity is illustrated graphically 
in Fig. 3. 

The information and the uncertainty about the density or specific 
gravity p to be expected for a meteoroid of mass m grams can be 
represented by 


0.108 +y 4 + 1.37 y -1.37 

p = 10 4 b m 

(95) 

y = -1.08 
7 4 

(93) 

y 5 = -0.079 

(97) 

cr = 0. 073 

(94) 



(r =0. 075 

(99) 


This mass dependence of density is illustrated graphically in 
Fig. 4. 

A spherical metallic shell (or empty tank) of density p t and of 
Brinell Hardness Number which, when in a near-earth orbit, must 
just resist being punctured when hit by a meteoroid of mass m grams, 
has a necessary wall thickness p centimeters, which is sufficiently 
accurately represented as having an approximately normally distributed 
common logarithm with the following mean and standard deviation. 

E [logjopl = 1.06 - 0. 29 log 10 P t - 25 l°gioH t + 0. 338 log 10 m 

(134) 
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^logiop ~ {°- 053 + [°- 204 + 0. 056 log 10 (p t /H t ) + 0. 0046 log 10 m] 2 + 

+ [o. 093 + 0. 086 log 10 p t + 0 . 0068 log 10 m ] 2 + 

+ [o.OZl + 0.0151og 10 mj Z J 2 (135) 


Two examples of metals which might be used in such manne 

(Pt»Ht) = (2.80, 135): hard aluminum alloy 

= (7.42, 310): hard stainless steel, 

and these are the basis for the results which are illustrated 
in Figs. 5 and 6. 


r are: 

(136) 

(137) 
graphically 


The number of punctures per second per square meter of exposed 
hemispherical area of a vehicle in a near-earth orbit can be considered 
to have an approximately normally distributed common logarithm with 
mean and standard deviation: 

E [ lo gio4>] =-9.13-3. 52 logjo p - 1.01 log 10 p t - 0.881 log 10 H t 

(139) 

^logio’t’ ~ fl • 8 6 + (0. 1 61 logic p - 0.0956 + 0. 0460 log 10 p^. + 

+ 0.0401 logic H t ) 2 + (0 . 0474 logic p + 0. 666 + 

+ 0. Ill logic p t - 0. 0859 logic H t ) 2 + (-0. 0709 logic p - 

- 0.252 - 0.323 log 10 p t - 0 . 01 77 logic H t ) 2 ] 2 (140) 

These results for the two metals described above are illustrated 
graphically in Figs. 7 and 8. 


The number of centimeters p of equivalent free wall thickness, for 
a vehicle of hemispherical area A square meters exposed in a near- 
earth orbit for t seconds, necessary to give a no-puncture probability 
R, can be considered to have an approximately normally distributed 
common logarithm with mean and standard deviation given by: 

E [logic p] = 0.284 [log^ At - logic (-log e R)]- 0. 288 logic p t - 

- 0.25 logic Hf- - 2.37 (1551 
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logio p 


■{°- 


150 + 


0. 00382 (logio At - log 10 (-log e R)) + 

+ 0.0556 log 10 (pt/Ht) + 0 . 1 9 8 J + 

+ Jo. 0057 (log l0 At - logio (-log e R)) + 0. 086logi 0 p t + 

+ 0 . 0 1 9 6j 2 + Jo. 0129 (log l0 At - log 10 (-log e R)) - 

YY (156) 


- 0.145 


These results are illustrated graphically, in rigs. 9 through 11 for the 
hard aluminum alloy described in Eq. 136, and in Figs. 12 through 14 
for the hard stainless steel described in Eq. 137. 


The variance has been resolved into nine independent components, 
the numerical relative magnitudes of which are presented in Table 1. 

Some further sources of uncertainty must be considered in 
practice, e. g. , effects of contents and composite— wall structures. 


SECTION I. INTRODUCTION 


A. SCOPE 

Since about two years ago several excellent reports [e.g. , 3, 

6, 12, 13, 17, 19, 25, 29, 33, 34, 35, and 4 l] have been directed 
specifically toward meteoroid hazards. None of these has dealt with 
the problem using operations research and systems analysis techniques. 
These techniques have been effective against similar difficulties in 
military and other situations [e.g., 7 and 1 1] . 

In the five stages described by Way, et.al. [40^ , in which 
reports may be most appropriately designated, this report is best 
described by "In the third stage, review articles appear. These put 
together data and theories and often come up with a working model, 
a handy description, ..." 

The report is intended to be more than an illustration, and 
more than a mathematical model. Single- shell-of -metal is the basis 
of the analysis. Of course, there is considerable interest in composite 


5 



bumpers and, e.g. , in multiple- sheet structures [e.g. , 33] ; but it is 
hoped that the effectiveness of other materials and structures can be 
considered in terms of the effectiveness of some thickness of single- 
shell- of -metal . Similarly special considerations must be given for 
vehicles which: (1) have attitude control or (2) do not remain in a near 

earth orbit or (3) do not remain near the ecliptic plane. 

B. METHOD 

The value that a physical parameter (e.g., length, velocity, 
etc.) will have under specified circumstances will be uncertain when: 
(1) the parameter (or the physical process of which the parameter is 
indicative) may be capricious (e.g. , the parameter may be a random 
statistical variable) , (2) the available information may not be either 
sufficiently firm or direct, or (3) the latter is compounded with the 
former - as in the present problem - the resulting uncertainty being 
treated as if it were due to randomness. 
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SECTION II. MATHEMATICAL PRINCIPLES 


A. FORMAL BASIS 

Let x x , . . . , x n be any continuous and statistically independent 
variables with probability density functions f(x x ), . . . , f(x n ). Then, 
because of statistical independence, the joint probability density function 
can be expr e s s ed [e . g . , 24] by 

f(x lf .... x n ) =f(x 1 )... f(x n ) (1) 


One of the most pertinent of the mathematical statistical 
concepts is the mean or expected value of a function g(x x , . . . , x n ) of 
the chance variables x x , . . . , x n ; i. e. , 


30 


oo 


E |g (xj , . . . , x n ) 


/ ... J g(x lf . . . , X n )f(x x x n)dx x . . • dx n . 

( 2 ) 


- co 


By Eqs . 1 and 2, as the x x , . . . , x n are statistically independent, 


E ax "f . . . a n x n 


K 


= a, E x, + . . . + a n E 


[ x n 


- a l X 1 + • * • + a n x n 


(3) 


and 


E jx x . . . x n j - E [x x ] . . . E [x n ] xj • • • x n 
where 3 l \ , ...» a-n are constants. 


(4) 


A. convenient measure of the dispersion of a continuous chance 
variable x is its standard deviation c r x , the square root of the second 
moment with respect to the mean [e. g. , 24] ; i. e. , 


oo _i 

0" x = [/(X- x) 2 f(x) dxj ( 3 * 5 1 

-00 

By Eqs. 2 and 5, one can define the standard deviation of a random 
variable x as the square root of the expected value of the square of the 
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deviation of the variable x from its expected value x. Therefore, by 
similarly defining ctq, the standard deviation of a function G (x* , . . . , x ) 
of the statistically independent chance variables x 2 , . . . , x n , and using 1 
Eqs. 1 and 2 


00 00 

•c =[/•••/ (c - E[G]) Z f(x 1 )...f(x n )dx,...dx n ]? (6) 

-00 -00 


But Eq. 6 is usually not sufficiently convenient for finding 
numerical results. The formula, previously argued pragmatically by 
the author [8], is the same as that which has, for the special case of 
normal variables, been rigorously derived by Scarborough [36]; i.e. , 



A more general proof of Eq. 7, for variables not necessarily 
normal, can be given. When the information about the variables is 
given by discrete data (e. g. , with the values x u , . . . , x llM - for x, ), then 
for Eqs . 1 and 5, 


Xl ~ ^ (Xu + ' * ' + x iN> (8) 

**1 = {b [ (XU " Xl )Z + ' • • + (x iN ~ x * >']}! (9) 


But by taking the total differential of G(x x , . . . , x n ). 


dG 


9G 

0X X 


dxt + 


+ — dx n 
9x n n 


(10) 


Equation 10 is approximately true when all of the differentials are 
replaced by sufficiently small increments. Assume that the 
deviations from the mean in Eqs. 8 and 9 are sufficiently small incre- 
ments.* Then corresponding to the N sets of data for the n variables 
one has, by substituting into Eq. 10, the N equations 

Gi " G " 9^ (Xli " Xl) + * ' • + < x ni ■%) * (11) 


where i - 1, . . . , N. Therefore, by squaring both sides of Eq. 11, 
adding the N resulting corresponding equations and dividing by N, 

* A more accurate general formula than Eq. 7 would be obtained by 
including non-linear terms in Eq. 11. 
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N 


*.Z < G i ' G) 


Nitl 




N 

^ " ^j>< x Ui -%) 

where 


(12) 


Z ( x ji " XjMxki - Xk) = 

1N i = 1 

1 N i N i N i N __ 

= F I XjiXki-Xj -ft 1 Xki - Xk- Z Xji +- Z XjXk • 

i = 1 i-l i-l i-l (13) 


The last term on the right side of Eq. 13 is identically xj x^; and by 
Eq. 8, the two middle terms are each equal to -5cj5q^; while the first 
term is the discrete case expression for the expected value of the 
product of two independent variables xj and x^ and, therefore by Eq. 4, 
is also equal to xj x^. Therefore, Eq. 13 vanishes and, by Eq. 9, that 
which is left of Eq. 12, after taking the square root of both sides, is 
identically Eq. 7. 


After the probability density function for a chance variable has 
been decided, then other interesting results can be derived; e. g. , 
consider a normal variable y with mean y and standard deviation c r, 

(1) the probability density function is [e. g. , 24] 


y 


f (y) = (2ircr y 2 ) 


i -(2a v 2 )-My - y) 


— \2 


(14) 


(2) there is an even chance that a randomly selected sample of the 
variable y will have a value within the (probable error) interval 


Y - 0. 6745cr y £y < y + 0. 6745<r y (15) 

and (3) the expected value of e a ^, where a is a constant, by Eqs. 2 and 
14 is 
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E[e ay ]= Je ay f ty) dy = 
-00 


= /< Z^)-K'^y y ' [ y ‘ <V +ai7 y 2 )] ! +ay + ja 2 o\ 


ay + 1 a 2 <r 2 
e y 


' dy 
(16) 


More generally than by Eq. 15, when x is any continuous 
random variable, not necessarily normally distributed, with distribu- 
tion function f(x), then the probability C (confidence) that 


x * x 6 

is related to x^ and Xj^, the lowest possible value of x, by 

x 6 

C = J f (x) dx 
X L 


(17) 


(18) 


where 6 is an index by which a particular value of x is established. 
When x is a differentiable function of a normally distributed random 
variable y; i. e. , 


x = x(y) ( 19 ) 

then, symbolically or otherwise, by solving Eq. 19 explicitly for y as 
a function of x: 


y = y( x )- (20) 

By applying the technique of Eqs. 2 and 14 to Eq. 19 , one can find 
E [x] , the expected value of x, to which there corresponds a value of y 
by substituting E [x] for x in Eq. 20. Then it should be anticipated that 
the most appropriate general functional relation between 6 and y^ is 

y 5 = y( E M) + 6cr y (21) 

where: (1) Eq. 15 is a special case, and (2) the x§ upper limit in 
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Eqs. 17 and 18 is found by substituting the right side of Eq. 21 for y in 
Eq. 19- By introductory mathematical statistics [e.g. , 24], the f(x) 
probability density function for x in Eq. 18 is found by: (1) substituting 

the right side of Eq. 20 for y in the Eq. 14 expression for the f(y) 
probability density function for y, and (2) multiplying by the absolute 
value of the derivative of the right side of Eq. 20 with respect to x. 
Therefore Eq. 18 becomes 
x (y(E [x]) + 6cr y ) 

C = J (2no- y 2 ) 

X L 

This Eq. 22 establishes the general functional relation between 
confidence C and index 6, the Eq. 15 results being only for a special 
case; i. e. , the special case being x = y in Eq. 19- 


i e -( 2 o y Z) ~ l y( x ) - y] 2 r.d_ 


L dx 


y (x) 


dx. 


( 22 ) 


Another special case with which one necessarily must be 
concerned here is where, except for a constant coefficient, x in Eq. 19 
takes the following functional form: 

x = (23) 


where (3 is a constant. Then, by Eqs. 2, 14, and 23, 

p(y + i|3<r ! ) 

E[x] = e y 

and, by Eqs. 20 and 23, 

y = (p‘ l ) log e x. 

Also, by Eqs. 24 and 25, Eq. 21 becomes: 

y 6 = y + iP°- y 2 + 6<r y . 

Therefore, by Eqs. 21, 23, 25, and 26, Eq. 22 becomes: 

rt / — i 1 r» Z t C \ 


(24) 


(25) 


( 26 ) 


P(y + iP°y + 60 


C = (2irp 2 cr y 2 )~ 2 


7 


( x ) e 


~(2cr y 2 )' 1 [(p)- 1 log e x-yj ; 


dx . 


(27) 
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So, by Eq. Z7, a definite value of C follows from a definite value of 6, 
the indicated integral having become a definite integral, and x having 
become a dummy variable - which is related to a more appropriate 
dummy variable t by: 

P(y + t) 

x - e 7 

with which transformation one can write Eq. Z7 in the following 
convenient form: 



-00 


where the values of the indicated integrals are available from tabulations 
[?•§•> 24] of normal areas. But with the tables of normal areas one can 
show that Eqs . 23, 26, and 28 lead to the following convenient result: 

y y- 0.6745cr y y + 0.6745cr 

(C, 10 ) = (0.25, 10 X ), (0.50, 10 ), (0.75, 10 y ). 

(29) 

B. OPERATIONAL METHODS 

What has been done in an uncertain situation? Considerable 
information (more or less accurate and more or less pertinent) has 
been painstakingly obtained over many years and at great expense; but 
it is not amenable to easy and certain application to the meteoroid 
hazard problem. Some say there is probably a considerable hazard 
ahead; while others say possible but not very likely. A design decision 
is upcoming. If the hazard is real, then plans for counteracting it must 
be made promptly. But if the hazard is illusory, then counteractions 
will have been an unnecessary bother and expense. 

Efforts to gain the much needed new information are being 
appropriately encouraged; e.g. , Whipple [42j is said to indicate that 
80% of the needed meteor information may be obtained by terrestrial 
observation. The results of such efforts will facilitate and improve 
design decisions which eventually must be made for later spacecraft; 
but the results will not be available for the present purpose. Reviewing, 
calculating, and discussing are the only efforts for which there is yet 
sufficient time; and as a by-product of these efforts one desires also to 
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learn how most efficiently to expend the forthcoming information collecting 
efforts . 


In other words, the meteoroid hazard to space vehicles is a 
typical example of the problems for which operations research methods 
have been developed. Of course, Ehricke [l9] "The [meteoroid flux] 
numbers, therefore, vary by several orders of magnitude. On this basis, 
a reliable estimate regarding collision probabilities and erosion effects is 
obviously impossible" and Rodriguez [35] "The most certain aspect of 
the hazard to space vehicles from meteoroids is the great uncertainty 
underlying any attempt at a quantitative assessment of the problem" are 
right, but more pertinent is Singer's [38] comment "It is generally 
conceded that measurements of meteoric particles are among the most 
difficult and uncertain of satellite experiments. But this makes it all 
the more incumbent on the investigators to report their findings m the 
greatest detail possible". 

C. PROCEDURAL CONVENIENCE 

One is justified in going ahead and trying to use to some 
advantage whatever information can be found regardless of how meager. 
For example, the following formula derived in [10] relates confidence C 
that reliability is not less than R when n t tests have given n f failures 
when n f is a very small but not necessarily vanishing part of n t : 

R=[e- n£ (l - C)] <nt ' ^ < 30 > 

where 0 ^ C ^ 1 when rif = 0, and where £ - C ^ 1 when nj + 0. Surely 
with no failures from two tests it is helpful to think that, based only on 
demonstrated performance, there is an even chance that the reliability 
is in the interval 0 . 50 — R — 0. 87; or in some other interval determined 
by hypothetically increasing n t in Eq. 30 [e.g. , 1 1] to make allowance 
for information other than test performance. More specifically, what 
this has to do with procedural convenience is that Poisson statistics 
were assumed in the derivation of Eq. 30; whereas it is not necessarily 
strictly true that the failure process in the above example is Poisson. 

But it was tacitly assumed as a reasonable approximation. 


So a procedural convenience is an approximation which is made 
as an appropriate expedient, without which the analysis would be unaccep- 
tably difficult or even impossible. Other examples which could have been 
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given for illustration are, e.g. : ( 1 ) assuming variables to be normal 

when in fact they may be only very nearly normal, ( 2 ) assuming chance 
variables to be statistically independent when in fact they may have 
some small correlation, and (3) assuming a curve to be piece-wise 
constant by representing it by a histogram. But as Davison and Winslow 
[l3j have stated "The common as sumption[in space vehicle hazard 
considerations] is that Poisson statistics describe the [meteoroid] 
population", and with appropriate reservations, "... the assumption of 
Poisson statistics as descriptive of the population is not completely 
unreasonable. " 


There appears to be some further need for interpreting the 
use of Poisson statistics as mentioned in several of the reference 
reports [l-42]. Consider an event which must happen some time and 
which can happen any time. Then, as Toralballa [39] has shown, the 
probability density distribution for the occurrence of the event is 

t 

-/ m(t) dt 

f(t) = m(t) e ° /31 


where m(t) is a function of time which in life statistics is called 
"specific mortality" and in engineering statistics is sometimes called 
"failure-rate-of-survivors". Then the product m(t) dt is the conditional 
probability that if the event has not already occurred by time t then it 
will occur in the interval between t and t + dt. Then, by the author's 
derivation shown in a previous report [ 9 ], the probability that the event 
will not have occurred by time t is 


R = 1 


t 



f(t)dt = e 


0 




(32) 


where (in meteoroid technology) <j>A is flux <j> times exposed area A. Any 
procedural convenience to be gained in the determination of R in Eq. 32 
will depend essentially on the convenience of evaluating the integral of 
m(t); but, as Toralballa [39] has shown, the specific mortality function 
corresponding to the Poisson distribution function is the constant 
reciprocal of the mean-time -to-occurrence for the repeated event; i.e., 
m is the average rate of occurrence of the event, or in meteoroid 
technology, f lux <|> times exposed area A, and Eq. 32 becomes 


R = 


-4>At 

e 


(33) 
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SECTION III. INTERPRETATION OF THE AVAILABLE 
INFORMATION 


A. METEOROID FLUX 

1. Temporal Dependence of Flux . As McCoy [29] states 
"historically, celestial debris has been classified as either 'shower' 
(swarms and streams) or 'sporadic 1 . Any meteor that could not be 
identified with a stream or swarm was classified as sporadic. For 
meteors in the visual range, several investigators have shown a 
preponderance (60-95%) to be sporadic; Lovell has estimated 80-90% 
sporadic for all meteor sizes. 11 

Davison and Winslow [l Z] state that shower meteors make 
up about 20-30% of the meteors sighted; and they cite several references 
to the effect that "The peak activity during these showers may be 4 to 
5 times the sporadic meteor rate, but on extremely rare occasions 
much greater. " 

Ellyett and Keay [18] have shown some data which may 
indicate a significant variation in hourly meteor rate throughout the 
day and in daily rate throughout the year "... expected as a result 
of the tilt of the earth's axis, causing a hemispherical autumnal 
maximum and vernal minimum. " Also Dubin [15] in an analysis of 
meteoritic impact data on artificial satellite Explorer I (1958a) says 
that, after discarding 16 of the first 88 hits, "A diurnal dependence, 
resulting by the earth's heliocentric velocity, is evidenced, in that 
nearly 90% of the 72 hits shown . . . occurred on the dawn side of the 
earth between the hours of midnight and twelve noon. " Hibbs [22 and 
23] has subsequently given a different interpretation of the data 
(Section III. A. 3), but failed to convince Dubin [16]- 

Bel'kovich [4] in reporting analytical and experimental 
investigation of the Poisson distribution hypothesis for "... the chance 
recording of a meteor, at a radar station or during visual observations 
..." concluded: "The distiibution of the number of meteors in time 

follows the Poisson law and is independent of the size of the observa- 
tional time interval and of the presence of showers. " In his same 
report Bel'kovich [4] also said: "The distribution of the number of 

recorded meteors in time will not be stationary since observational 


15 


conditions are continuously changing owing to the earth's daily rotation 
and motion along its orbit. " Perhaps what Bel'kovich means by the 
above two statements is that the rate will change but that one can predict 
neither when the change will occur nor what the rate will be after it has 
changed. 


But in Eq. 33 different values of <f> may be forecast for 
different time intervals, as by histogram representation. One can 
then see that a non-Poisson situation can be approximated as piece - 
wise Poisson, for cj> is piece-wise constant. Then, instead of Eq. 33, 

R = e" A [ 4 ’ l(tl ' W + • • • +4fa (*n - *n - i >] (34 

for t = t n - t 0 . 


It should be emphasized, however, that Eqs. 33 and 34 
give the same result for the same exposure At whenever <j> in Eq. 33 
is the time average of the <h. • • • ,4> n in E q- 34. It does not make any 
difference whether or not the $ in Eq. 33 and the 4> x , .... 4> n in Eq. 34 
are further resolved into: (1) components due to sporadic meteoroids, 

and (2) components due to shower meteoroids - unless it should be 
found that the two populations should be separately described with 
respect to mass, density, velocity, etc. 

2 - Directional Dependence of Flux . Dubin [j.5] writes that 
Because the earth is shielding the satellite, the impact rate for the 
same mass component of micrometeorites at one astronomical unit is 
nearly twice this [detected] value. . . Whipple [4l] writes with 
appropriate caution: "A correction for distance above the earth's 

surface is rather difficult to apply unless the precise orbital character- 
istics of the incoming meteoric bodies are known. Perhaps the inclusion 
of the factor i is adequate for the ordinary problem near the earth's 
surface. At extreme distances, greater than 10 4 km, a complete 
recalculation is needed with very careful attention to the orbital char- 
acteristics of meteoroids in space and of the space vehicle in question. 
Even after correction for the factor j, the number striking the vehicle 
will probably fall off with Increasing distance from the earth. " 

According to Siedentopf [37]: "The particles entering the 

earth's atmosphere [meteor observations] have very eccentric orbits, 
whereas most particles that are optically effective move in nearly 
circular orbits under the action of the Poynting- Roberts on effect. The 
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observable meteors generally have radii greater than 100p, whereas 
the brightness of the zodiacal light and the Fraunhofer corona comes 
mainly from the particles with radii between 1 and 100p. " Also 
Ehricke [19] says: "According to Takaknbo the radius of the dust 

particles most effective in the zodiacal light and the F corona is about 

20p. " 


McCoy [29] reports: "These sporadic particles enter the 

earth's atmosphere from all directions, within approximately ±25° of 
the ecliptic, and at irregular intervals. It has been observed through 
analysis of radiant points that the ratio of meteors approaching the 
earth from behind to those approaching from ahead is close to 30:1. 


This is because most meteors have direct orbits like the earth rather 
than retrograde motion. However, the true ratio of direct to retrograde 
orbit meteors is approximately 50:1. The difference in ratios is 
explained by the fact that the earth overtakes some of the slower meteors. 
Also, because of the earth's orbital velocity more meteoric encounters 
(~2:1) occur on the leading or forward face, as can be readily substan- 
tiated by observing the increased number of visual displays after 
midnight than before. M 


Davison and Winslow [13] show a graph (which they credit 
to Lovell) for the distribution of visual sporadic meteor radiants in 
ecliptic latitude which supports the estimate that approximately 50% 
are inclined within ±15° of the ecliptic. And in two Other graphs (which 
they credit to Hawkins [51]) they [12 and 13] show: (1) "Polar diagram 

drawn in the plane of the earth's orbit which shows the apparent number 
of meteor radiants detected per unit angle, per unit time" and (2) "Polar 
diagram drawn in the plane of the earth's orbit which shows the number 
of meteors per unit angle which cross earth's orbit per unit time. 

Those two graphs, reproduced from Hawkins' [5l] paper, are included 
in this report as Figs. 15 and 16. They show clearly the predominant 
proper motion of the meteoroids and the tendency for the earth to over- 
take the meteoroids. Davison and Winslow [12] rightly remark, It 
could be concluded, therefore, that the preferred orientation of a 
[spacecraft's exposed] surface is parallel to the plane of the ecliptic 
with the major axis parallel to the apex of the earth's way. " 


If the spacecraft is to have attitude control, then part of 
its area, Aj , will be so located and oriented that the meteoroid flux 
with respect to it is <J>j , and similarly for A 2 with c |> 2 . etc. ; so that, 
instead of Eq. 33, it would follow from Eq. 32 that 


17 



-t(4>iA 1 +... +4’n A n). 


(35) 


R = e 

3 - Spatial Dependence of Flux . Nazarova [3l] concluded: "If 

one compares the results obtained from experiments with Sputnik III 
and the three Soviet space rockets, it will be seen that the density of 
meteoric matter in the vicinity of the earth is not constant. It varies 
in time and space. " 

After some analysis of the micrometeoroid-impact sonic 
data from Sputnik HI (1958 6 2 ), Nazarova [32] concluded "... that the 
number of impacts varies with changes in the position of the satellite 
in its orbit but not with changes in altitude. " And subsequently, in 
interpreting the micrometeoroid-impact sonic data from Explorer VIII 
(I960 £), McCracken, et.al. [30], similarly concluded ".. . that avail- 
able direct measurements show no discernible evidence of an appreciable 
geocentric concentration of interplanetary dust particles." 

But that meteoroid flux should be accentuated in the vicinity 
of a planet, aside from the shadowing effect discussed in Section III. A. 2. , 
should be expected because it can be shown analytically that, for hyper- 
bolic trajectories of meteoroids passing planets, the partial derivative 
of radial distance of closest approach with respect to the projected (i. e. , 
f ree) miss distance is both positive and less than unity. In 
regard to this phenomenon, the author has gone no further than to this 
qualitative verification, but Beard [3] writes: "The overall effect of the 

dust blanket at the surface of the earth will be to increase the concentra- 
tion of dust, calculated in the absence of the earth's gravitational field, 
by a factor possibly as large as 6000, depending on the eccentricity of 
the dust's orbits. (The factor [Beard says] is more likely about a few 
hundred.) . . .A third of this blanket is concentrated within one earth 
radius of the earth's surface ... A simple approximate calculation of 
the dust orbit in the gravitational field of the sun and a planet reveals 
that the concentration of the dust varies inversely with the three-halves 
power of its distance from the center of the planet. " But one suspects 
that the radial dependence of meteoroid flux must also be functionally 
related to the mass of the particular planet. 


That meteoroid flux in the vicinity of a planet may be 
radial-distance dependent is further suggested by capture processes 
other than direct collision. In 1958 Baker [2] reported: "Preliminary 

calculations show that only about 0. 2% of the total number of the porous 
stony meteorites which strike the earth will result in natural satellites 
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[which]. . . would originate from ’near-miss' meteoritic trajectories 
that only graze the atmosphere of the earth, the meteorites being slowed 
sufficiently to enter onto a geocentric elliptical orbit. . .It is noted that 
there are other [further] mechanisms for capture of natural satellites 
involving the attraction of the moon, the Poynting-Robertson effect, and 
so on ...” By including data which had been earlier discarded by 
Dubin [15], as discussed in Section IH.A.l. , Hibbs [22] gave this inter- 
pretation to the micrometeoroid-impact data from artificial satellite 
Explorer I: "The distribution in longitude relative to the satellite - 

earth- sun angle corresponds to an altitude distribution and apparently 
contains no information that is not better shown in this latter distribu- 
tion. With suitable analysis, the altitude distribution yields information 
on the velocity of the particles relative to the center of the earth. The 
conclusion is that the average particle measured by Explorer I was in 
a closed orbit around the earth rather than on an impact trajectory from 
a great distance to the surface of the earth. " As each captured particle 
would repeat many orbits, supposedly the 0.2% captures, quoted above 
from Baker [2] , would be sufficient to support a considerable steady- 
state population of captured particles. These findings by Hibbs are 
further supported by those later reported by Laevastu and Mellis [27] : 

11 The estimates of the rate of fall of cosmic spherules to the earth, 
based upon the counts from sediments, depend, among other things, on 
the exact determination of the rate of sedimentation. As the value 
calculated from satellite data is considerably higher, it gives added 
support to the conclusions of Hibbs. . . " 

Aside from the above hypothesized increase in meteoroid 
flux in the vicinities of the planets, there are further observational 
evidences that meteoroid flux increases with decreasing distance from 
the ecliptic and from the sun. Beard [3] says: ’’The dust diffracts 

sunlight at small angles which may be observed in the solar corona 
during solar eclipses. Sunlight is also reflected by the dust at large 
angles which may be observed as . . . zodiacal light. ,! This interpre- 
tation of solar corona and zodiacal light data is further supported by 
the in-fall theory of cosmology; e.g. , (1) Beard [3] says: ’’The 

relativistic interaction of these particles with sunlight, as shown by 
Robertson and Wyatt and Whipple, causes them to spiral slowly into 
the sun with a radial velocity inversely proportional to the solar distance 
. . . One further factor that determines the radial distribution of inter- 
planetary dust is the gravitational attraction of the planets. . . which 
causes the dust to approach the plane of the ecliptic with a drift velocity 
that is inversely proportional to the square root of the solar distance. 
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Inclusion of this effect results in a dust concentration in the ecliptic 
plane that is inversely proportional to the three-halves power of the 
solar distance, " and ( 2 ) Best [5] says: "Supposing the solar system 

encountered such a cloud [of interstellar dust at a relative velocity 
corresponding to the 20 km/sec solar proper velocity], the particles 
would all have hyperbolic orbits with respect to the sun. However, the 
orbits do not remain hyperbolic of fixed eccentricity due to the braking 
force of solar radiation, which for initially elliptical orbits leads to the 
Poynting -Robertson effect, and the eccentricity will decrease near the 
sun . . . Note that [the criterion for capture by retardation by solar 
radiation] is identical with [the criterion for the Lyttleton cometary 
accretion process] and thus if either capture process is involved the 
capture radius will be about 5 A. U. , which is just inside the orbit of 
Jupiter. " 


4. Mass Dependence of Flux . So as not to have to be 
concerned further in this report with the directional aspects of 
meteoroid flux, one can assume that the term is used in the same 
sense as Whipple [41] used it: the meteoroid flux encountered by a 

sphere in a near-earth orbit is the number per unit time hitting the 
sphere divided by the area of the effectively exposed hemisphere. And, 
as Rodriguez [35] suggests: "An equation of the form 

F > = (36) 

is used by many investigators. F^, is the [flux] of particles having 
mass m or greater, and a x and pj are empirical constants." 

The selection of a solution pair of values (c^ , Pj ) is not 
decided in the (F^ , m) domain by most authors, excepting Laevastu 
and Mellis [27], It is decided in the (logjoF > , logjo m) domain; i.e. , 
from a family of equi -probable contours of possible solution points 
which, through Eqs . 16 through 29, are the consequences of the 
following underlying functional representation 


logic^ = Px log 10 m + Yl ( 37 ) 

where Pi is a constant and Yi is an approximately normally distributed 
random variable with mean yj and standard deviation c Ty . This Eq. 37 
and such other functional representations as will be used in this report 
to represent uncertain quantities and/or relations are intended as exped- 
ients for the propagation of confidence through other functions of this 
and other uncertain parameters. 
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Functional relations between the flux of meteoroids, 
with masses equal to or greater than m, to be encountered by a space- 
craft in orbit near the earth, or relations between log 10 F> and log 10 m, 
have been inferred primarily from data obtained by careful quantitative 
measurements of physical parameters involved in the following four 
phenomena: (1) the interaction of meteoroids with the atmosphere of 

the earth as studied by visual, photographic, and radar methods and 
discussed (e.g.) by Whipple [41], (2) the accumulation of meteoroid 
debris on the earth as studied by chemical, physical, and statistical 
analyses of the sediments and discussed (e.g.) by Laevastu and 
Mellis [27], (3) the disturbance by micrometeoroids hitting instru- 
mented artificial satellites as studied by microphones and discussed 
(e.g.) by Dubin [l4 and 15], and (4) interaction between micrometeoroid 
and electromagnetic radiation as inferred by physical optical peculiar- 
ities of solar corona and zodiacal light and discussed (e.g.) by Beard [3]. 
The combined interpretation of these four separate categories of infor- 
mation is further facilitated (or confused, depending on one's point of 
view) by theoretical considerations of the paths of meteoroids , of 
given mass and cross-sectional area, moving under the combined 
influences of solar radiation and solar and planetary combined gravita- 
tional fields - as discussed (e.g.) by Siedentopf [37], Best [5], and 
Beard [3]. 


The various published interpretations of the information 
described above differ considerably, and in other ways also suggest 
that the results must yet be uncertain to a considerable extent. When 
various interpretations are considered as points in a graph, logjo F^ 
versus log l0 m, then results are missing for 10‘ 7 gm < m < 10' 4 gm. 
Except for the degradation of optical surfaces, thin films, paints, etc. , 
this is the entire region of interest - mass m possibly large enough 
for penetration of a space vehicle (or vital component) and flux F^ 
possibly too high to be ignored. Also not a single one of the various 
phenomena or underlying principles of interpretation gives points on 
both sides of the interval of interest. 

Concerning visual, photographic and radar tracking of 
meteors, Nazarova [32] says: "At present [i960] . . . these methods 

are only useful for registering particles with masses of 10 4 gm and 
more." For this side of the region of interest, Whipple [41] suggests: 
". . . an increase in the number of meteoroids inversely as the mass of 
the particle. " This is equivalent to assuming a unit value of -Pi in 
Eq. 37. 


21 


( 38 ) 


By taking the antilogarithm of Eq. 37 one finds 

F - Pi ( lo 8e 10) y! 

■t - m e 1 

which (except for the constant factor) is the same functional form as 
Eq. 23 with 


(3 = log e 10 = 2. 303 . (39) 

Therefore Eqs . 23 through 29 are applicable here; and Whipple's [4l] 
tabulated data imply that 

?l = " 12-20 - *f x (40) 

where ^ is a component which would correct any bias in m 0 in Eq. 42. 

About the micrometeoroid side of the region of interest, 
Bjork [6] says: ". . . At the present time [May I960], rocket and 

satellite experiments are inherently limited to measuring the impacts 
of meteoroids having masses of about 10~ 8 gm or less." By McCracken, 
Alexander, and Dubin's [30] more recent interpretation of Explorer VIII 
results, Eq. 37 is approximately true for lO' 10 < m <10’ 6 gm when 


Pi = -1.70 


and 


Y x = -17.0 


B J ork f6] recommends a relation, to fit both the satellite 
and radar data regions, which similarly implies that 

Pi = -10/9 


and 


?i - - 12 • 

Nysmith and Summers [33] gave a summary of results from 
ten different satellite determinations of (3 X and y x in Eq. 37 with mass 
sensitivities in the interval 10' 10 < m £ 2 x 10‘ 8 gm and with an average 
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value of -8. 70 for the common logarithm of the mass sensitivity limit. 
They preferred a unit value for -(3 X , for comparison with meteor data, 
with which, by equally weighting their tabulated results, 

y = -11. 28 
y l 


<r v = 0. 60 . 
y i 

Beard [3] finds, from studies of solar corona and zodiacal 
light, that flux depends on radius to roughly the 3.5 power for small 
particles; implying that: 

(3 X = -3.5/3 = -1.17 . 

The establishment of a mathematical model which bridges 
the intermediate interval of extrapolation is not to be accomplished 
without considerable reservation. Laevastu and Mellis [27] say. In 
making this extrapolation [from satellite data on small particles to 
larger particles], account is taken of the fact that essentially all cosmic 
particles > 250p in diameter are intercepted by the planet Jupiter and 
therefore do not reach the earth. " But since the volume of a sphere 
with 250 (jl diameter is 8.2xl0' 6 cubic centimeters, the mass must be 
within the interval 4x 10~ 7 ^ m ^ 3x 10“ 5 gm for density (Section III. C.) 
within the interval 0.05 ^ p — 3.6. Also, about the small particles 
affecting solar corona and zodiacal light, Siedentopf [37] says: M A 

comparison of these results with the data obtained from meteor obser- 
vations is difficult. The particles entering the earth's atmosphere 
have very eccentric orbits, whereas most particles that are optically 
effective move in nearly circular orbits under the action of the 
Poynting-Robertson effect. The observable meteors generally have 
radii greater than 100p whereas the brightness of the zodiacal light 
and the Fraunhofer corona comes mainly from the particles with radii 
between 1 and IOOjjl. So the information gained by optical methods and 
by meteor observations covers different aspects of the interplanetary 
matter that practically do not overlap." Interestingly enough, these 
apprehensions support the notion that for the high-flux region, say 
F > > I0- 7 particles per square meter per second, the results of 
McCracken, Alexander, and Dubin [30] may be appropriate; whereas 
Whipple's [41] results may be appropriate for lower flux values. This 
would be equivalent to assuming that there are no meteoroids with 
masses in the interval 
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10'( I7 ° - 7 °>A 7 = 1 . 3 x 10’ 6 < m < io-l n -88-7.oo) - 


- 1. 3 x 10" 5 gm (41) 

because they have been captured by the major planets. It appears to be 
more likely that the available information on either side of the region 
of extrapolation is not sufficiently firm to warrant interpretations such 
as Eq . 41 . 


Now consider the accuracy of the data for the larger 
meteoroids. As the basis for his analysis, Whipple [4l] says, ". z 
meteor of visual magnitude zero is determined to have a mass of the 
order of 25 gm, " Bjork [6] says that Whipple . . has calculated mo 
lies between 1 and 30 gm with 25 gm the preferred value, " 


In the physical theory of meteors, as presented by Opik 
f44] , Kallmann [45] , Levin [46], and Jacchia [47], the expressions for 
meteoroid mass m, as related to visual magnitude, are always found 
by taking the antilogarithm of a linear expression for log 10 m, for which 
the slope and intercept are both somewhat random and problematical 
and, the present state of the information being such as it is, might 
just as well be assumed to be normally distributed. This is part of the 
basis for the contention that yi can be considered to be approximately 
normally distributed; but also it shows that values of m 0 bear the Eq, 29 
relation to confidence. Then the 25 gm value for m 0 is a result of 
choice rather than of skewness; the value otherwise implicated being 


2 (l°gi o 30 + logic 1) 

m 0 =10 


5.48 . 


(42) 


discussing IvlcCrosky's interpretation of his experiments 
using simulated meteors to estimate luminous efficiency, Davison and 
Winslow [12] stated McCrosky indicated that luminous efficiency from 
3 to 100 times as large as that used by Whipple was not precluded by 
the experiment. The common logarithm of such a correction factor 
would be the further bias which would have to be algebraically 
subtracted fromyjin Eq. 37; because Whipple [41] says: "It is 

assumed that the mass decreases by a factor of 10 04 per magnitude 
step ..." while a unit increase in visual magnitude corresponds also to 
a luminous intensity decrease by a factor of 10 04 . Therefore, the bias 
in y t corresponding to a luminous efficiency correction factor of 25/m 0 


24 



a Yj = logio (25 /m 0 ) 


(43) 


The question of what value to use for m 0 is related to the 
km /sec meteoroid velocity v a by Opik [44] as follows: 

log 10 m 0 = 10.97 - 1-7 log 10 (10 5 v a ) - 


This Eq. 43 implies that m 0 is 1 . 02 or 0. 70 gm depending on whether 
one uses 28 km/sec for v a , as Whipple [41] does for all visual magni- 
tudes up to 7, or 35 km /sec, as will be shown in Section III. B. to be 
more nearly deducible from mass versus velocity relations. 


When all of the above reasons are considered, the author 
prefers to use the Eq. 42 value for mo; then by Eqs. 40, 42, and 43, 


y =-12.86. 

7 1 


(44) 


Also by the reasons already discussed in this section, it 
follows that when -8. 70 is the mean of the common logarithm of the 
mass sensitivity limits for the satellite data, as discussed above, then 
the median mass for the detected micrometeoroids will be approximately 


io -8.70 + (-(3 1 )' 1 log 10 2 


gm = m s 


(45) 


at which value for m one has also implied that 

logi 0 F > | = 8.70 + (P 1 )- , log 10 2 - 11.28 = 


= -2. 58 + 0. 301 (f3j )-' . (46) 

But by Eqs. 37 and 44 one has also implied that, at the same value of 
m, 


logio F > | = Pj[-8. 70 - 0. 301 (pi)' 1 ] - 12. 86 = 

m s 

= -8. 70Pi - 13. 16 . (47) 

The simultaneous solution of Eqs. 46 and 47 is 
Px = -1. 19 (48) 
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which agrees very well with the -1.17 value, mentioned above, based 
on Beard’s [3] estimate. 

Nothing has been said about the confidence value which 
should be as sociated with the above mentioned delimitations : 1 <m 0 < 30 gm, 
but presumably it is intended that they are each separated from the 50% 
confidence value for m 0 by not less than one probable error of m 0 ; i.e. 

cr y 1 = 0.50 (log 10 30) (0 . 6745)- 1 = 1.10. (49) 

In deciding whether this Eq. 49 estimate for cr^ is too much in excess 
of the value 0. 60 mentioned above, one must consider further the various 
uncertainties besetting the determination of the masses of micro- 
meteoroids based on data from instrumented satellites and from other 
sources . 


The mass dependence of micrometeoroid flux for m^lO' 7 gm, 
as inferred from microphone impulses from instrumented satellites, 
also has considerable possibility for error. In reporting Explorer I 
results, Dubin [l 5] says: "It was also assumed that since the micro- 

meteorite detector on 1958a* was sensitive to the momenta of the impacts 
to velocities as high as 4 km/sec, the momentum transfer remained 
proportional to the first power of the velocity of impact at meteor velo- 
cities. If this is not basically correct, the error may become large at 
the high end of the velocity distribution ... It was further assumed that 
all impulses measured were meteoritic impacts.” Elsewhere, Dubin [14] 
says. The velocities of meteoroid impacts upon a satellite may vary 
between 10 km/sec and 70 km/sec. . . " and says further that, according 
to computations by Bjork, at 70 km/sec the vehicle experiences a change 
in momentum which is ”... two and one-half times the initial particle 
momentum.” More recently McCracken, Alexander, and Dubin [30] 
report that the calibration has been extended through speeds up to 
approximately 10 km/sec and that: "An average velocity of 30 km/sec, 

i dative to a satellite has been assigned to the dust particles in order that 
the data may be used to determine the average mass distribution. . .” 

And, in inferring the mass of micrometeoroids from the impulses from 
microphones in Sputnik III, Nazarova [31 and 32] used a similar 
relationship - . • while assuming that the mean velocity of the particle 
is equal to 40 km/sec. . . ” But Whipple [41] says: "Undoubtedly the 

velocity falls off for smaller meteoroids as we deal more and more 
with particles whose orbital eccentricities and dimensions have been 
reduced by physical effects. . . " and for m < 10~ 7 gm he assumed an 
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average velocity of 15 km/sec. Also Opik [43] infers that a velocity 
distribution with mean velocity and standard deviation equal to 14.5 
and 2.4 km/sec respectively implies a frequency distribution of particle 
diameters which is in good agreement with the distribution of diameters 
found by Laevastu and Mellis [43] for ’’extraterrestrial material in 
deep-sea deposits. M Opik [43] concluded that: M A low value of and a 

low dispersion in the velocities is definitely indicated, corresponding 
to cosmic dust circling the sun in direct orbits of low eccentricity. ,f 
It does seem possible, therefore, that any bias due to a proportionately 
larger transfer of momentum above the calibration velocities may 
compensate for any bias due to assuming higher average velocity. But 
with so much room for doubt, it seems appropriate to accept the Eq. 49 

criterion for cry . 

y l 

The author prefers to consider that, throughout the 
interval 10" 10 < m < 10 1 gm, meteoroid flux F > and mass m are related 
by Eqs. 29, 37, 44, 48, and 49 as illustrated graphically in Fig. 1 
based on the Eq. 29 criteria. 

B. METEOROID VELOCITY: MASS DEPENDENCE 

1. Relative to the Earth's Atmosphere. Whipple [41] says: 

"A velocity of 28 km/sec is average for photographic meteors. 
Undoubtedly the velocity falls off for smaller meteoroids. . . The 
velocity at the edge of the earth's atmosphere cannot be less than 
1 1 km /sec because of the earth 1 s attraction ... A mean value of 1 5 km / sec 
has been arbitrarily chosen for the smaller meteoroids and an arbitrary 
gradation of velocity with magnitude adopted.” Thus, in his table, 
Whipple [4l] indicated a velocity of 28 km /sec for all visual magnitudes 
less than or equal to 7, 15 km/sec for all visual magnitudes greater 
than or equal to 20, and otherwise a 1 km/sec decrease in velocity for 
each unit increase in visual magnitude. Also, in Section III. A. 4. it 
was shown that Whipple [41] assumed that the mass decreases by a 
factor of 10°* 4 per magnitude step, and therefore he implies a linear 
relation between velocity and log 10 m, which is apparently the basis for 
the illustration given by Nysmith and Summers [33] . In his tabulation, 
Whipple [41] further shows 3 . 96 x 1 0“ 10 < m ^ 25 . 0 gm for 27< magnitude =£ 0 . 
But in the discussion leading up to Eq. 44, the decision was that the 
masses which Whipple assumed should be reduced by the factor 10~ 0,66 ; 
the above mass interval being thereby changed to 10" 10 ' 06 ^ m^l0°* 74 gm. 
But if Whipple [41] had used, between magnitude 27 and 0, the same 
linear relation between magnitude and velocity which he assumed 
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between magnitudes 20 and 7, then the velocity would have been reduced 
nearly to the circular orbital velocity at magnitude 27, in agreement 
with Hibb’s [22 and 23] analysis of micrometeorite data from instru- 
mented satellites, and at magnitude 0 the 35 km/sec velocity would 
agree very nearly with the 35. 5 km/sec average for the data for larger 
meteors attributed by Davison and Winslow [l2 and 13] to Hawkins and 
Southworth. This would imply the following relation between the 
velocity, v a (km/sec) and mass, m(gm): 

v a - (35 - 8)(logio 10° 74 - log 10 lO -10 06 )' 1 log 10 m + 

+[8 - (35 - 8)(log 10 10° 74 - log 10 1 0- 10 06 )-> (log 10 lO' 10 06 )] 
= 2.50 logxo m + 33. 2 (50) 

for 10‘ 10 ^ m £ 10 1 gm. 

Dubin [l 5] refers to the cosmic dust detection limit of 
optical and radio measurements as n . . . approximately visual magni- 
tude 10 ... " One wonders what the corresponding limit would be for 
the determination of meteoroid velocity at the edge of the earth’s 
atmosphere -- evidently the magnitude limit would not be quite as high 
as it is just for detection. Whipple [41] says only: ’’Radio meteor 

studies may provide measures of the [velocity] gradation in the near 
future, perhaps to the twelfth magnitude. 11 

By adjusting Whipple’s [41] tabulated data for meteoroid 
masses, as discussed above, the relation between visual magnitude M 
and mass m is 


in 0. 74 - 0. 40M 
m = 10 


(51) 


One can consider the average velocity of all meteoroids between the 
relatively infrequently occurring magnitude 0 and some higher magnitude 
M u 5 10 as a function of M u . By Eqs . 37 and 50, 

(3. (0. 74 - 0. 40M) + y, 

= 10 1 (52) 


Then by Eqs. 48 and 51, the number of meteors with M < M u is greater 
than the number with M < M u - 1 by the following factor: 


-0. 406, 
10 hl 


3.00 . 
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Therefore the average velocity for meteors of magnitude * M u will be 
very nearly the same as the average velocity for those of a magnitude 
in the vicinity of M u . On that basis one would expect that for Whipple's 
[41] data M u is not appreciably higher than 7. No corresponding 
estimate is available for the above-mentioned data which gives velocity 
mean and standard deviation 35.5 and 13.3 km/sec respectively, but, 
assuming the Eq. 50 corresponding value of m, one has: 


(m, v a , <r v ) 


(10 


0.92 


gm, 


35.5 km/sec, 13. 3 km/sec). 

(53) 


Opik's [43] theoretical distribution of cosmic dust 
velocities and diameters below 300p gives velocity mean and standard 
deviation 14. 5 and 2.40 km/sec respectively, and an average value 
of -6.99 for logio (m/p). By taking the value of log 10 m which, by Eq. 5 , 
gives the 14.5 km/sec average velocity, one finds 


(log 10 m, log 10 P) = (-7.48, -0. 49 = logio 0.32) 


(54) 


i. e. , p = 0. 32 for small meteoroids. This agrees well with other 
results; e.g., Davison and Winslow [12] say: "Based on Beard’s 

estimates [based on solar F -corona and zodiacal light] it would appear 
that the material density of these particles [lp <diameter < 300pJ must 
be equal to or less than 0.3 gram per cubic centimeter. " So one can 
agree that there is sufficient basis for another set of values: 


(m, 


v ■ 


) = (10 -7 ' 48 gm, 14.5km/sec, 2. 40 km/sec) . 


(55) 


But the following Eq. 5 6 functional relation 
logio v = P 2 logio m + y 


(56) 


where p 2 is a constant and y 2 is an approximately normally distributed 
random variable with mean y 2 and standard deviation oy^ will more 
appropriately represent the Eqs. 53 and 55 information, than will one 
like Eq. 50, because: (1) ^y 2 decreases as v a decreases, and (2) values 

of v 3 must not be too severely depressed as m decreases. 

CJL 

By applying Eq. 7 to Eq. 56 and using the Eqs. 53 and 
55 information, 
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a y 2 " <r io gl0 v. 


= {<r v a y/v a ) log 10 e 


(13. 3/35. 5) log! o e = 0 . 1 63 ^ 

I independent estimates 
(2.40/14. 5)log 10 e = 0.072 ^ mean: 


= | (0. 163 + 0. 072) =0.12 

B Y Eqs. 29 and 56, it follows also from the information of Eos 53 
and 55 that ^ 

log 10 35. 5 = p 2 log 10 1 0° 92 + y 


(57) 


lo gio 1^* 5 — p 2 log 10 1 0‘ 


. 7 . 48 , - 


Yz 


Then by simultaneous solution of Eqs. 58 and 59, 
P 2 = 0.046 


(58) 

(59) 

(60) 
(61) 

The results from Eqs. 57, 60. and 61 for <r p 2 , and y ; respectively 
when used m Eq. 56 with the criteria from Eq. 29 for contours of 
equi- probable values of v a versus m, seem to represent very well the 
available information throughout the interval 10- 10 < m < 10 1 em as 
illustrated graphically in Fig. 2. s 


y 2 =1.51 


2 . 


Relative to a Vehicle in a Near-Earth Orbit. As the 


- — ** X 1^0.1 - L11 uru u, t 

earth orbits around the sun, the gravitational field of the earth 
deflects toward the earth the meteoroids which otherwise, in their 
e iptical orbits around the sun, would merely pass close by the earth. 

The distribution of closing velocities between the 
meteoroids and a space vehicle in orbit near the earth will depend on 
the distribution of angular deviations of the local vertical from the 
tangents to the meteoroid trajectories at the point of impact. It is 
assumed: (1) that the result will be approximately the same if a two- 

body meteoroid-earth model is substituted for the more nearly correct 
three-body meteoroid-sun-earth model, and (2) that, just as meteoroid 
ve oci y with respect to the earth's atmosphere can be related to mass 
by Eq 56, closing velocity v c with respect to a space vehicle in orbit 
near the earth can be similarly related by 


l°glov c - p 3 l°g 10 m + y 


(62) 
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where p 3 is a constant and y 3 is an approximately normally distributed 
random variable with mean y and standard deviation cm representing 

D J 3 

uncertainty of the relation physically and according to available 
information. 


The necessary functional relations for hyperbolic 
trajectories are given in introductory treatments of celestial mechanics 
[e.g., 19]. The hyperbolic velocity excess v^ is 


v oo = < v a 2 - 2 ^ m e 


a) 


(63) 


where m e is the mass of the earth, 


gravitation, and r 


y is the universal constant of 
is the radial distance from the center of the earth 


to the zone implicit in meteoroid velocity data discussed above. The 
velocity v r at radial distance r is, by Eq. 63 


v r = 


v a 2 - 2ym 


(64) 


Because a spacecraft is orbiting at geocentric radial distance r, it is 

convenient to consider that a geocentric sphere of radius r is making 

a tunnel of radius Di through a swarm of meteoroids, where. (1) the 

2 it 

trajectory of any meteoroid originally approaching along the surface 
of the tunnel with the closing velocity v^ will be tangent to the sphere r, 
and (Z) the trajectory of any meteoroid originally approaching with the 

same velocity but within the tunnel and displaced only D x ^ ^ rom 

1 . ^ 

the tunnel axis will cut the sphere r at angle xx < with respect to 
the local position vector (i.e. , x x is the zenith angle). Because of the 
conservation of angular momentum, the dot product of the meteoroid 
velocity vector and geocentric radius vector is invariant with respect 
to the position of any particular meteoroid along its trajectory; i.e. , 


v D v = V' r sin x, 
oo m 1 


(65) 


Therefore, by Eqs. 63 through 65, 


sinx 1 = (l/r)D 


*1 


r a v a 2 - 


, r a v a 2 “ 2 \ m e( 1 " r a/ r U 


= /Dj 


Iir. 


( 66 ) 


Then by Eq. 66, the tunnel radius is 
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Di 

IV 


1 ~ 2Ym e (l - r a /r)/r 
1 - 2Ym e /r a v a 2 


a v a 


1 

2 


( 67 ) 


Although the above tunnel concept was convenient in the derivation of 
Eqs. 66 and 67, it should be emphasized that the results are valid also 
when meteoroids approach the earth from all directions as they do in 
practice; e. g. , by Eq. 66 one sees that half of the meteors should 
approach the atmosphere with angles deviating less than 45 degrees 
from the zenith. 


When a spacecraft is orbiting in a nearly circular orbit 
with geocentric coordinate r, then its velocity v s is essentially 
horizontal with 


v s = (^ m e /r) 2 ■ (68) 

Then, when one assumes that both the meteoroid and space vehicle are 
moving in the same plane, it follows from Eqs. 64, 66, and 68 that 
the closing velocity v c is 

v c = [( v s + v r s ^ n x l ) 2 + (v r cos x x ) 2 J 2 =: 

= [ v s 2 + v / + 2v r v s sinx i] J (69) 

Because the relative area of a narrow concentric ring in the cross- 
section of the tunnel is equal to the product of its relative circumference 
and its differential relative width, and because one must admit equal 
probability for positive and negative values of Xj , it follows from Eq. 66 
that the probability distribution function for xj is 

f ( x i ) = | (2tt sin Xl )-£- (sinxj. )/2it | = \ |sin2 Xl | . (70) 

By Eqs . 2 and 70 

x : =0. 

By Eqs. 5 and 70 

°- Xl = [(^ 2 - 4 )/ 8] 1 

By Eqs. 2 through 4, 69, and 70, the expected value of v c 2 is 


(71) 

(72) 
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(73) 


E K 2 ] 


+ E[v r 2 ] . 


But because the closing velocity v c is increased from v r by only a 
relatively small increment due to v s , Eq. 73 can be appropriately 
approximated and, with Eqs. 64 and 68, simplified by 


v. 


(v s 2 + v r 2 ) 2 = V a 2 + (ym e /r)(3 - 2r/r a ) 


(74) 


By Eq. 7, by taking the standard deviation of v c in Eq. 69, and then by 
applying the relations from Eqs. 64, 68, 71, 72, and 74, 


r 

I 


9 v c 9v r 


cr 


9v^ 3vn v a 


+ 


9v c 

9x x 


\ 2 t 

%) J 


~ Rva *v) Z + ( v r v s ff x, 

v r L a 


/?= 

= (v a (Ty^) 2 + (v r v s [ v s v r j 

= [ V a 2(tr v a + <r x 1 2 ^ m e/ r) ' Z (y™ e / r)Z (r/: 
v a 2 + (\m e /r)(3 - 2r/r a )j_ 2 


- l)cr 


1 

2 12 


X 1 


(75) 


Levin [46] says that 6740 km is the radius of the atmospheric layer in 
which meteors are observed. Then for a space vehicle in orbit near 
the earth, say at (6378 + 500) km, 

11 1 
(ym e /r) 2 = (ym e /r a ) 2 = 7. 9 1 ( 6378 / 6878) 2 = 7.62km/sec . 

(76) 


With Eqs. 72 and 76, Eqs. 74 and 75 become 

v c = [v a 2 + (7 - 62) 2 ] 2 


V r 


■{ 


V 


+ (7. 62)' 


MF 




(77) 

(78) 


These Eqs. 77 and 78 results, together with Eq. 62 and the Eqs. 3 and 
7 operations, lead to the following Eqs. 79 and 80 relations: 


p 3 log 10 m + y 3 




(79) 
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°y 3 = a iogio V C = (v/ v c) lQ gio e = 

= v a (logioe) [cr Va 2 + (7. 6Z) 2 (tt 2 - 4)/8j 2 £v a 2 + (7. 6Z) 2 j- 1 . 

(80) 

The following Eqs. 81 and 8Z are found by substituting the information 


from Eqs. 53 and 55 into Eq. 79: 

p 3 log I0 10°- 92 +7 3 = log 10 [(35.5) 2 4- (7.6Z) 2 ]* (81) 

P 3 log 10 10" 7 - 48 +y 3 = log 10 [(14. 5) 2 + (7 . 6Z) z l 2 (8Z) 

By solving Eqs. 81 and 8Z simultaneously, one finds the following 
numerical values for P 3 and y • 

P 3 = 0. 041 ( 83 ) 

y 3 = 1 - 52 - (84) 


The following Eqs. 85 and 86 are found by substituting the information 
from Eqs. 53 and 55 into Eq. 80: 

(v a , <rv a> <r y J = (35.5, 13.3, 0. 173), (14. 5, Z.40, 0.168) 

(85) 

«r y = \ (0. 173 + 0. 168) = 0. 17 . (86) 

The results from Eqs. 83, 84, and 86 for p 3 , y 3 , and oy respectively, 

when used in Eq. 6Z together with the Eq. Z9 criteria for contours of 
probable values of v^ versus m, seem to represent very well the 
available information throughout the interval 10" 10 < m < 10 1 gm, as 
illustrated graphically in Fig. 3. 

C. METEOROID DENSITY: MASS DEPENDENCE 

Meteoroids differ not only in mass, but also those which have 
approximately the same mass may differ widely in composition and 
structure. They are, in the order of increasing abundance and 
decreasing puncturability , classified broadly as: (1) metallic, with 

density somewhat more or less than that of iron, say p = 7.8, (Z) stony, 
say with density somewhat more or less than p = 3.5, and (3) fluffy, 
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with density sometimes said to be somewhat more or less than p = 0.05 
Whipple [4l] says: M . . . the meteoric bodies, or meteoroids, mani- 

fested as ordinary meteors are extremely fragile and breakable. n So 
presumably none of the fluffy meteoroids of appreciable size would have 
sufficient integrity to survive through the atmosphere and be recovered 
from the earth's surface. Then one may further presume that it is 
only with respect to metallic and stony meteoroids that one may apply 
Holl's [26] comment: "From the analysis of nearly 1000 meteorites 

recovered from the earth's surface, it has been determined that 
dielectric particles are nearly ten times as numerous as metallic ones. 

Whipple [4l] says: "The evidence is extremely strong, 

although not quite conclusive, that the density of ordinary meteoroids 
is the order of 0.05 gm/cm 3 . This value depends, to a limited extent, 
upon physical arguments from photographic meteor data. . But one 
is somewhat at a loss for an interpretation of !! . . .is the order of. . 

This can more appropriately be decided after reviewing the discussion 
in Section III. A. 4. about the possibility that, as Davison and Winslow 
[12] attribute to McCrosky, the luminous efficiency for meteors may 
be from 3 to 100 times as large as that used by Whipple. Davison and 
Winslow [12] say further that: "An increase of 100 in the efficiency 

would have two ramifications. It would decrease the mass associated 
with a particular meteor sighting by a factor of 100 and increase the 
particle density for spherical particles by a factor of 10. " Apparently 
the functional relation between density p and mass m, which has to 
represent both physical fact and the extent of information about it, 
might be represented by: 


logloP = P 4 logi 0 m + y 4 (87) 

where (3 4 is a constant and y 4 is an approximately normally distributed 

random variable with mean y„ and standard deviations cr__ . Because 

y 4 y 4 

of the (25/5.48 = 4.56) meteor luminous efficiency correction factor 
(in the interval 3<4.56<100), the choice for which has already been 
supported in Section IH. A. 4. , Whipple's above estimate of the density 
of "ordinary meteoroids" must be transformed to 

1 

(v a , p) = [28 km/sec, 0. 05 (4. 56) 1 2 = 0. 107], (88) 


35 


But, by Eqs. 56, 60, and 61, the 28 km/sec velocity of ’’ordinary 
meteoroids” corresponds to the mass of the same "ordinary meteoroids” 
by: 


(v a , m) = (28km/sec, 


1(J (logio 28- 1.51)(0.046)- 1 


lO' 137 gm) . 
(89) 


Therefore, by Eqs. 88 and 89, the mass and density correspondence 
for ordinary meteors is: 

(m, p) = (1CT 1 * 37 gm, 0. 107 = IQ' 0 - 972 ) . (90) 


Because meteoroid density varies with the square root of the meteor 
luminous efficiency factor, and mass varies directly with it, and if 
one will assume in the absence of more definitive information that, 
at any specified value of mass, the standard deviation of log 10 p bears 
approximately the same relation to the standard deviation of log 10 
as it would if all of the variances in them were due only to an uncertainty 
in the meteor luminous efficiency factor, then, by Eqs. 7, 37, 44, 48, 

49, and 87, 


°Y 4 ' °log 10 P ~IM a log 10 m-|P4/Pi|cr logio - |p 4 /Pi| oy - 

= | ^/PilN = |1. 10/ -1. 1 9 1 1 p 4 | = 0. 924 |p 4 |. (91) 

Estimates for values of density p to correspond with mass m for the 
micrometeoroids are not so firmly established as for ordinary meteoroids. 
Bjork [6] says: "It is reasonable to expect that the smaller meteoroids 

will not have the lacy, porous structure needed to give a specific gravity 
of 0.05, but that their density will increase as one goes down the mass 
scale, eventually approaching the [2.8] specific gravity of stone, the 
basic component, as very small sizes are reached. " For the density of 
particles of 7p radius, Alexander, McCracken, and LaGow [1] as sumed 
1.4xl0‘ 9 gm, saying: "The value lgmcm -3 for mass density is some- 

what higher than the value 0. 05 gm cm* 3 used by Whipple. The value used 
by Whipple is probably much more appropriate to the photographic 
meteors, for which it was determined, than for the dust particles in the 
direct measurements range of sizes." But Whipple [4l] cautions: "We 

have no sound basis upon which to assume densities for smaller than 
ordinary meteor particles, except for the effect of light pressure or 
other physical effects. . . " However Opik's [43] analysis, as discussed 
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in Section III. B. 1. , leads to Eq. 54, which the author prefers to accept 
also for the purpose of the present report. Then as Eq. 87 is that of a 
line with two points given by Eqs. 54 and 90, 


p 4 = -0.079 (92) 

y = -1.08. (93) 

y 4 

Then also, by Eqs. 91 and 92, 


cr = 0.073. 

y 4 


(94) 


These Eqs. 92 through 94 results for p 4 , y 4> and oy respectively 

when used in Eq. 87 together with the Eq. 29 criteria for contours of 
equi-probable values of p versus m, seem very well to represent the 
avialable information when the mass m is in the vicinity of that 
( 10 ~ 1 , 37 gm) for M ordinary meteoroids”; and the 50% confidence contour 
seems appropriate over the entire interval 10" 10 < m < 10 1 gm . But the 
25 and 75 % confidence contours must diverge more for small values of 
m because of the greater uncertainty. Therefore, Eq. 87 will be 
replaced by the following more flexible model: 


logio P = y 5 (logiom - (3 5 ) + y 4 + (3 4 (3 5 (95) 

where: ( 1 ) ( 3 4 , y 4 , and cry are as before, with the values given in 

Eqs. 92 through 94 respectively, (2) P 5 is the constant common 
logarithm of the Eq. 90 mass of the "ordinary meteoroids" 

p 5 = - 1 . 37 (96) 


and ( 3 ) y^ is an approximately normally distributed random variable, 
statistically independent of y 4 , with standard deviation cry and mean 
y 5 . The latter is, by Eqs. 87, 92, and 95, 

= p 4 = -0. 079 - (97) 

Then by Eqs. 29 and 95 through 97, the 50% confidence contour of p 
versus m is invariant with respect to replacing Eq. 87 with Eq. 95. 
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The c r y standard deviation will be so chosen that the above 

7 5 

estimate by Alexander, McCracken, and LaGow [i]. which can be 
represented by 

(p, m) = ( 1 , 1.4x1 0 -9 - 10~ 8 ' 85 gm) , ( 98 ) 

is on the 75% confidence contour; i.e. , on the contour of values which 
have a 25% chance of being less than p. Therefore, with Eqs. 7, 29, 
and 95 through 98, 

o- y = 0 . 075 . ( 99 ) 

5 

By applying the Eqs. 7 and 29 criteria and the Eqs. 92 through 
94, 96, 97, and 99 numerical values to Eq. 95 one gets the contours of 
equi-probable values of p versus m which are illustrated graphically in 
Fig. 4. The results seem to be satisfactory throughout the interval 
10 - 1 ° < m < io 1 gm. 

D. METEOROID DAMAGE 

1. Nature and Function of Material Versus Effects. The 
meteoroid hazard to space vehicles, or the damage to be expected from 
meteors, is put in appropriate perspective by the following statement 
by Rinehart [34] about the concept of quality of failure: "In any partic- 

ular target, the failure will usually be a complex of many qualities of 
failure, although frequently a single quality predominates. The basic 
problems to consider are what qualities of failure prevail in the situa- 
at hand, which ones are of interest, and whether each quality of failure 
is an energy-absorbing process, a momentum -absorbing process, or a 
combination of both. And lastly, what is the quantitative relationship 
between the extent of failure of a particular quality and the energy and 
momentum available to cause the failure? ... a few overall qualities of 
failure are the following: perforation or puncture, volume of crater, 
volume of failed region, scabbing, spallation, and amount of abrasion. . . ” 

The quality of failure which will be emphasized in this 
report is failure by puncture of an empty metallic single shell. This 
has been necessary to keep the problem from getting too far afield. 

Should it be expected that the same shell would be more easily punc- 
tured when filled with a liquid? Or is it just the other way around? 

And then, of course, one would want to know what temperature, pressure, 
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shock wave, etc. , the liquid is subjected to in the temporal vicinity, 
what chemical reaction follows, and is there an explosive vaporization 
which further rends the structure? But such comprehensive considera- 
tions are beyond the scope of this report. 

2. Crater Volume in Thick Targets Versus Energy and 
Momentum for Meteoroids at Normal Incidence . Rinehart [34] illustrates 
his contention that the volume of the crater produced in a target material 
by an impacting meteoroid is a linear combination of the kinetic energy 
(| mv c 2 ) and momentum (mv c ) of the meteoroid with respect to the 
target. He says that: ". . .if a failure results from application of an 

impulse under which the material dislodges easily and offers the inertia 
of its own mass as a resistance to motion, then the process is a 
momentum transfer (a flicking away of material, so to speak). On the 
other hand, if the body steadily continues to resist application of the 
force, the process is energy-absorbing (pushing of material against a 
force). Perforation of a thin plate by a projectile, or penetration into 
a laminated structure such as wood, are momentum-absorbing 
processes. An energy-absorbing process is the formation of a deep 
crater in steel by the impact of a heavy projectile. In most real 
materials, the failure will be a combination of the two. " 

Both Eichelberger [17] and Beard [3] imply that, for any 
given target and meteoroid material and structure, crater volume is 
essentially proportional to kinetic energy ( j mv 2 ) independently of 
momentum (mv). Eichelberger [17], says: "The [empirical] results 

support very strongly the conclusion from fundamental considerations 
that cavitation plays the dominant role in crater formation in ductile 
materials and explains the linear relationship between volume and 
energy." Beard [3], says: "It is most probable that energy considera- 

tions of evaporation, rather than momentum effects, dominate the 
surface interaction of the micrometeorites with a satellite. " 

3. Crater Depth in Thick Targets Versus Energy , 

Momentum, and Density of Meteoroids . Typically the various empirical 
and theoretical formulas for meteoroid crater depth [e. g. , 48] are the 
product of the cube root of crater volume and the following three factors: 
(1) the crater-shape factor, (2) the plate -thickness factor, and (3) the 
angle- of- incidence factor. Thus, crater depth p 0 is proportional 
to the 1/3 power of the mass m if crater volume is a linear combination 
of energy ( \ mv c 2 ) and momentum (mv c ) . Since Bjork [6] indicates that 
p is proportional to (mv c )*^ he is therefore presumably implying that 
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cia.tci v olutTie is proportional to momentum (mv) independently of 
energy mv 2 ) . 


An example of a penetration law that does not infer that 
crater volume is a linear combination of energy and momentum, is 
the law that can be deduced from data which is graphically presented 
by Hoenig and Ritter [25]; i. e. , that the logarithm of penetration 
depth is linearly related to the logarithm of crater diameter in the 
same interval of kinetic energy over which they indicate with another 
graph that crater diameter is proportional to the square root of x 
kinetic energy. In other words, crater depth p Q is proportional to m 2 v 
for the experiments performed by Partridge at the University of Utah 
with wax pellets where impact velocity exceeds the velocity of sound 
in the target material. 


But to what must one say that the crater depth is 
proportional, when so many factors are available from which to choose? 
Interestingly enough the geometric mean of Hoeing and Ritter's [25] 
factor m 1/2 v c , and Bjork's [6] factor, (mv c ) 1/3 is (m 5/4 v c 2 ) 1/3 , which is 
very close to Eichelberger 1 s factor, (mv c 2 ) 1/3 . However, the author 
feels that neither momentum nor kinetic energy should be inconsequen- 
tial. He prefers, in the absence of more convincing contrary 
information, to consider that crater depth, at normal incidence and 
for given target and meteoroid material and shape, is represented by 


P 


0 



( 100 ) 


where: (1) (3^ is a constant which will be assumed to have the same 

value , 


^6 = 1/3 ( 101 ) 

as when crater volume is a linear combination of momentum and 
enei gy, and (Z) y^ is an approximately normally distributed random 
variable, which represents both randomness and uncertainty due to 
insufficient information. The algebraic mean for y^ is considered to 
be the exponent of the geometric mean of the factors, v c 1/3 and v 2/3 , 
which would represent proportionality of crater volume to momentum 
and energy, respectively, i.e., 

y 6 = ■ (102) 
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The standard deviation for y fe is considered to be one third of the 
difference between the above mentioned exponents 2/3 and 1/3, i.e. , 

1 / 3 < y^< 2/3 at approximately 87% confidence: 

<r = 1/9 < 103 > 

>6 

Herrmann and Jones [48] say: "Data on cratering has been 

reported by Summers for copper projectiles impacting copper targets at 
7, 000 and 11, 000 ft/sec [2. 13 and 3. 35 km/ sec], and by Kmeke for steel 
discs impacting lead targets at 16,400 ft/sec [5.00 km/ sec]. Both 
experiments noted that the data for oblique impact compared very well 
with that for normal impact, if penetration versus the normal component 
of velocity is plotted." This result can be represented by 


v 




v c COS X 2 


(104) 


where x 2 is the angle of incidence relative to the normal to the target 
surface and is therefore a random variable. Presumably this relation, 
Eq. 104, will also be sufficiently appropriate for meteoroids hitting 
other metal targets. By Eq. 104 the angle- of -incidence factor by which 

the right side of Eq. 100 must be multiplied is (cosx/ 6 , and therefore 
Eq. 100 is replaced by 


■ m 


1/3 


(v c cos X 2 




(105) 


To establish the statistical definition of x 2 in Eq. 105, 
consider a meteoroid incident on a sphere and an axis parallel to the 
path of the meteoroid but containing the center of the sphere. The 
plane normal to the axis and containing the center of the sphere divides 
the sphere into two hemispheres, and the one which is hit is orthogo- 
nally projected onto its base plane. If x 2 is the angle of incidence of 
the meteoroid relative to the normal to the surface of the hemisphere, 
then all meteoroids which are parallel to the axis and which will have 
angles in the interval between x 2 and x 2 + dx 2 will be projected onto a 
ring with radii in the interval between R sinx 2 and R sinx 2 + d(R sinx 2 ) 
The relative area of the differential ring is the probability density 
function for angular incidence x 2 , i.e., 


f(x 2 ) = |sin2x 2 | . (106) 

Since f(x 2 ) by Eq. 106 is symmetric with respect to itr , 
it follows that both the mean and median of x 2 are given by 
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x 2 = 4 if radians . 


( 107 ) 


Then by substUulmg th e Eqs. 106 and 107 expressions into Eq . 5, one 

finds that the standard deviation of x 2 is 


= (tt 2 /16 - |) 2 = 0.34 radians. 


(108) 


, Herrmann and Jones' [48] survey and analysis of published 

nlat empirical results for the crater depth in thick metal 

plates by hypervelocity projectiles are believed to provide a sufficient 
asis for establishing the crater shape factor. They conclude that, 
above a velocity transition region which depends on the target and 

-I* 6 Pr ° jectile stren § th d <>es not affect penetration, 
Lrto' , 7 ? projectiles , and: "If this is true, then the only 

actor to account for differences in penetration in a given target 
material by different projectile materials is the projectile density. " 

By an analysis of the empirical results published for many different 
target and projectiie material s> Herrmann and Jones [48] developed 
the following empirical non-dimensional penetration law 


P 0 /d = (0.36± 0.07)(p p /p t ) 2 /3( PpVc 2/H t )V3 


(109) 


where: (1) H t is the Brinell Hardness of the target, ( 2 ) p and p. are 

e projectile and target densities respectively, and ( 3 ) d is the 

Tl te + r ° f i hepr0jeCtile - The velocity exponent in Eq. 109 will agree 
with that m Eq. 100 when the last term in Eq. 109 is raised to the 
3y 6 /2 power, then Eq. 109 can be replaced by 


P 0 /d = 10 y7 (p p /p t )Z/3 (ptVc 2 /H t ) 


i 


( 110 ) 

where y ? is ; an approximately normally distributed random variable 

«:sr erainty in 1116 rela,ion and in 1116 about i.. 


10 y? = 0. 36 ± 0.07 = x 3 ± 0 -. 


x 3 


( 111 ) 


where one assumes from various comments that the numerically indi- 
cated uncertainty is standard deviation rather than probable error or 
mean deviation (mean absolute error). Also, Herrmann and Jones [48] 
have explained that they got the results, Eqs. 109 and 111 , by plotting 
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and fitting on log-log paper. So it is evident that x 3 in Eqs. 109 and 
111 is actually the antilogarithm of the more basic random variable y ? . 
Because it is subject to the normal law of error more directly in the 
fitting process than is x 3 , y 7 can more appropriately be considered 
as approximately normally distributed. Presumably y 7 can be considered 
to have mean y ? and standard deviation <Ty^ given by: 

Y v = logic 0. 36 = -0.44 (112) 

0 - = log 10 (0.36 + 0.07) - log 10 0.36 = 0. 08 (113) 

y 7 

But, when the units for p t are gm/cm 3 and those for v c are km/sec, 
then a further proportionality constant is necessary in the last term in 
Eqs. 109 and 110 because, by the definition in the Metals Handbook [50] 
the units for Brinell Hardness H-t are kilograms of force per square 
millimeter, i. e. , Eq. 1 1 0 must be replaced by 

p„/d = 10 y, (p p /p t ) 2/3 (kp t v c 2 /H t ) 2y ‘ (114) 

where 

(gm/cm 3 )( 10 5 cm /sec ) 2 - 102 ( 115 ) 

(980, 665 gm cm/sec z ) / (10 * 1 cm) 

Let both p and d in Eq. 114 be measured in centimeters. 
Then, assuming that meteoroid hazard is not essentially misrepre- 
sented by a spherical meteoroid of mass m grams. 

d = 2 ( 3 m/ 4 irpp ) l/3 . (116) 

Herrmann and Jones C48] make the following further 

comment about Eq. 109: "It might be noted that a slightly higher 

exponent in p^. might be expected to fit slightly better. However, it was 

decided to retain the advantages of a non-dimensional fit. Small 

changes (± 10 %) in exponents of the non-dimensional parameters did 

not significantly alter the mean deviation. " But because, by Eq. 102, 

y is 1 , the mean effect on p t of having adjusted the exponent of v c to 

y 6 is that the exponent of p t has been lowered. Also, the value to be 

used for the exponent of (pp/ p^) is uncertain, to a considerable extent, 

as is well illustrated by Bjork [ 6 ]. Therefore, an approximately 

normally distributed random variable exponent y g will be introduced 

with mean y , which is related to y, in that it restores the mean value 
8 ® 
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of the exponent, and with standard deviation cy indicative of a probable 
error which is 10% of the mean; i. e. , 8 


-y 8 + iy 6 = -2/3 + 1/3 

y 8 = iy 6 +1/3 = 7/12 = 0.583 (117) 

cr y 8 = o. 10(7/12) /0. 6745 = 0. 086 ( 118 ) 


Therefore, by replacing the exponent of the density ratio 
in Eq. 114 by y 8 , by substituting the Eq. 115 value for k in Eq. 114, by 
multiplying the resulting equation by Eq. 116, and by multiplying the 
side of the resulting equation by the angle-of— incidence factor 

y y 

(cosx 2 ) as in Eq. 105, one still gets the somewhat biased formula 
for the thick-target crater depth p Q due to a spherical meteoroid: 


10 y7 (102) 2y6 (2)(3/4 1 r) 1 / 3 P t 2Y6 ’ Y «h ’^P " 1/3 


m 


1/3 


(v c cos x 2 ) 




(119) 


when the following units are used: p Q , cm; p t and p gmcm' 3 ; H t , 
Brinell Hardness Number; m, gm; v c , km sec' 1 ; and x 2 , radians. 


Equation 119 is still biased because the slope of logjo p 
versus log 10 v c has been changed; but the intercept still ha« to be ° 
adjusted so that the combined effect of a change in slope and in inter- 
cept will be that some favored point is invariant. Eq. 11 9 . is based on 
Eq. 109, which is based on the upper segments of broken linear fits to 
data points in the log 10 (p/d) versus log 10 v c domain, for various 
projectile and target materials including hard aluminum alloys, and 
for log 10 v c typically in an interval of roughly log 10 3 km/sec < log 10 v c 
<l°glo 5 km /sec. Then take the favored point as 


(p p , p t , H t , v c ) = (2.80, 2.80, 135, 4) . (120) 

y 

Then the factor 10 9 by which one must multiply the right side of 
Eq. 119 to remove the bias is that which, at the Eq. 120 point, gives 
the following equality: 

y 1 

(P P /P t } 8p p l/3 {102 P t v c 2 /H t ) 2Y6 10 y9 = ( Pp /p t ) 1/3 (102v 2 /H t )»/ 3 . 

( 121 ) 
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y 

Then by Eqs . 120 and 121, the factor 10 9 is related to y 6 as follows: 

l/3 — ■ 

10 Y9 = [102 (4) 2 / 1 35] 2Y6 (2.80) ,/3 (122) 

Therefore, the above-mentioned bias is removed by multiplying the 
right side of Eq. 119 by the right side of Eq. 122, getting: 

p = 4. 01 (10) y7 (8. 44 p /H t ) 2V6 p t Ys p Y ® 1 m 1/3 (v c cos x z ) Y6 . 

° (123) 

It is of some interest to see how Eq. 123 compares with 
the following Eq. 124 for aluminum projectiles and thick aluminum 
targets reported by Bjork [6]: 

p = 1. 09 (mv_) 1/3 (124) 

0 c 

where the units are as in Eqs. 119 through 123. Assuming relatively 
pure aluminum metal with the following density and Brinell Hardness 
Number : 


(p t , H t ) = (2. 70, 23 <H t < 44), say (2. 70, 34) (125) 

then at normal incidence, with vanishing x 2 , one can equate the right 
sides of Eqs. 123 and 124, substitute the mean values 0.50, -0.44, 
and 0.583 for y , y , and y from Eqs. 102, 112, and 117 respectively 

and solve for the value of v c , i. e. , 

v c = 2.3 km/ sec . (126) 

Then the Eq. 123 formula gives values of thick-target crater depth 
larger than by the Eq. 124 formula by a factor (v c /2. 3) 1/6 , which is 

(6. 3/2.3) 1/6 =1.18 (127) 

at the 6. 3 km/sec value for v c used in the experiments on which the 
Eq. 124 formula is based. This Eq. 127 result is well within the 
Eq. Ill indicated experimental accuracy; and one does not know for 
sure if the Brinell Hardness Number of the aluminum was as low as 
is assumed in Eq. 125. The Eq. 127 factor vanishes when H^ is 
increased to 66, or to some value less than 66 if the density is also 
increased. One will therefore decide that Eq. 123 is satisfactory 
for thick targets. 
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4 - Thickness of a Just-Penetrable Shell Versus Mass , 
Density, Velocity, and Angle of Incidence of Meteoroids . The 
thickness p of a just-penetrable shell is related to the thick-target 
crater depth p Q by the target thickness factor 10 yi ° , i.e. , 

p/p Q = 10 10 (128) 


where y 10 is an approximately normally distributed random variable 
which represents both the randomness of the process and the uncertainty 
in the information about it. Bjork's comment is: "The calculations 

were made for thick targets, but enough information was obtained to 
deduce that if a projectile penetrates a depth p in a thick target, it 
will just penetrate a sheet of the same target material which is 1 . 5p 
thick. " Black [49] says: "To allow a 'bulge' although 'just not 

perforated', a skin gage of 1 . 5 times crater depth is generally assumed. 
(Note that this should be 2 - 3 if the results of Jaffe and Rittenhouse are 
used.)" Eichelberger's [1 7] "rule of thumb" is that: "... a meteoroid. . . 
will produce a hemispherical crater of volume t. . . If the thickness of 
the skin is less than (3t/2tt) 1/3 (or even if it is slightly greater), the 
skin will be perforated. " In other words, the factor is 2 ,/3 = 1 . 26 or 
slightly greater. Herrmann and Jones [48] illustrate experimental 
results (which they attribute to Kinard et.al. at NASA Langley) and an 
empirical formula indicating that, as shell thickness is decreased 
toward the value p, the value of p/p Q approaches: 


p/p -(1/1.3) 2 + 1 = 1.59 


(129) 


These Eq. 129 results are said to have been obtained with . . steel 
and aluminum projectiles into aluminum targets at impact velocities 
between 5, 000 and 13, 000 ft/sec" ... i.e., between 1 . 5 and 4 km/sec 
When C is the confidence that this Eq. 128 factor is equal to or less 
than the stated value, then it seems appropriate that the above- 
mentioned four estimates should be accepted as follows: 

(C, IQ 710 ) = (0.16, 1.26), (0.41, 1.50), (0.50, 1.59), 


(0.976, 2.5) 


(130) 


y 


10 


0.20 


(131) 
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oy =0.10. 
y ie 


( 132 ) 


SECTION IV. DESIGN AND OPERATIONAL PARAMETERS 


A. JUST-PENETRABLE METEOROID MASS VERSUS THICKNESS, 
DENSITY, AND HARDNESS OF THE FREE-WALL OF A 
SPACE VEHICLE IN A NEAR EARTH-ORBIT 


By taking the logarithm of the product of Eqs. 123 and 128 
and substituting the Eqs. 62 and 95 expressions for log 10 v c and log 10 P p 
one finds 

logio P = logio 4. 01 + y i0 + y 7 + iy 6 lo gl0 8. 44 + ly 6 i°gio (Pt /H t> - 


- y 8 logio P t + Ye lQ gi° < cos x 2 ) + y 3 y 6 + 

+ (y 8 - i/3)[y 4 + P 5 (P 4 -y s >] + [(> -v,)/3 + 

+ y 5 y 8 + P 3 y 6] logl ° m 


(133) 


With the values 0. 041, -0.079, -1.37, 1 j 1 _08, _ 

-0 079 i -0.44, 0.538, and 0.20 for (3 3 , (3 4 , P 5 > x 2 » Yz* Y 4’ Ys> Ye* 
y y , and y from Eqs. 83, 92, 96, 107, 84, 93, 97, 102, 112, 117, 
and 1 8 31 respectively, it follows by Eqs. 3 and 4 that the expected 
value of logio p in Eq. 133 is: 


E [logio p] = 1-06 - 0.29 logio P t - 0.25 logio H t + 0. 338 logio 


m . 
(134) 


And with the further values 0. 34, 0.17, 0.073, 0.075, 1/9, 0.08, 0.086, 
and 0. 10 for cr x , oy , oy^ <T y ^, oy^, oy^, oy^ and oy^ from Eqs. 108, 

86, 94, 99, 103 2 , 113, 118, and 132 respectively, it follows by Eq. 7 
that the standard deviation of logi 0 p in Eq. 133 is: 

p- p = -[o. 053 + [o. 204 + 0. 056 logio (P t ^ H P + °* 0° 4 61og IO m ] 2 + 

+ [o. 093 + 0. 086 logio P* + 0. 0068 logio mj 2 + 

- ..1 

+ [o. 021 + 0. 015 log 10 mj z | 2 . 


(135) 
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These Eqs. 133 through 135 results are illustrated graphically 
in Figs. 5 and 6 for the following two respective metals: 

(Pf = (2.80, 135): hard aluminum alloy 

= (7.42, 310): Harder stainless steel 

and for the assumption that log 10 p in Eq. 133 is sufficiently nearly 
normally distributed that the Eqs. 134 and 135 values for the mean and 
standard deviation can be used in the Eq. 29 formula relating confidence 
and values of the antilogarithm of an approximately normally distributed 
random variable. 

B. METEOROID PUNCTURE-FLUX VERSUS THICKNESS, 

DENSITY, AND HARDNESS OF THE FREE (EMPTY) WALL 

OP A SPACE VEHICLE IN A NEAR-EARTH ORBIT 

The penetration and/or puncture flux <j> is that value of F in 
Eq. 37 corresponding to a jus t-punctur able value of m. The following 
Eq. 138 formula for <J> is found by solving Eq. 133 explicitly for log 10 m 
and substituting the resulting expression into Eq. 37; i.e., 

logic 4> = + ^{logjo p - log 10 4.01 - y iQ - y? - iy 6 log 10 8.44 - 

- iy 6 iog X o (P t / H t ) +y 8 iog 10 p t - y 6 iog 10 (cosx 2 ) - 
" y 3 y 6 - (y 8 - 1/3 )[y4 + P 5 (P4-y5)]}[(i-y 5 )/3+y 5 y 8 + (3 3 y 6 ]->. 

(138) 

The random statistical variable 4>, the puncture-flux is found, 
by taking the antilogarithm of Eq. 138, to be a function of nine statis-’ 
tically independent random variables (y lf x 2 , y 3 through y 8 , and y ), 
all of which are approximately normally distributed except x 2 . If the 
right side of Eq. 138 were just a linear combination of these random 
variables and if x 2 were nearly normally distributed, then log 10 4> would 
also be a normal random variable, with mean given by using the means 
of the other random variables in Eq. 138, and Eq. 29 would be rigorously 
applicable for finding the contours of equi-probable values of c|> as a 
function of p. The difference here is that the right side of Eq. 138 
involves some products and ratios rather than just a linear combina- 
tion of the random variables. But it will be assumed to be a sufficiently 


(136) 

(137) 
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accurate approximation here that log l0 4> is normally distributed with a 
mean value given by using the means of the other random variables in 
Eq. 138. Then, with the values -1.19 and -12. 86 for and ^ from 
Eqs. 48 and 44 respectively, and with the other numerical values that 
were used in Eq. 133 to get Eq. 134, it follows from Eq. 138 that the 
mean or expected value of log 10 4> is: 

E flogxo*] = -9. 13 - 3.52 log 10 P - 1-01 log 10 P t - 0.881 log 10 H t . 
L J (139) 

And with the value 1. 10 for from Eq. 49, and with the other 

numerical values that were used to get Eq. 135 from Eq. 133, it 
follows by Eq. 7 that the standard deviation of logio^ i n Eq. 138 is: 



+ 0. 0401 log 10 H t ) 2 (147) 
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2 




(0. 0474 logjo p + 0.666 + 0. Ill logjo - 
- 0. 0859 log 10 H t ) 2 (148) 


( 9 0y^ 10 ^ j = (-0- 0709 log 10 P - 0. 252 
8 

- 0.0177 logjo H t )f 


0. 323 logjo P t - 


(149) 


These results, Eqs. 138 through 149, applied to each of the 
two metals described by Eqs. 136 and 137, lead to the following 
formulas, Eqs. 150 through 15 3, for the means and standard deviations 
of the common logarithms of the puncture -fluxes <|> as a function of 
shell thickness p: (1) for the hard aluminum alloy, 

E [l°glo4>] = -3.52 1og 10 p - 11.46 

°’logio<l> = [ 2 ‘ 34 + °‘ 146l °Sio P + 0. 033 (logjo p) 2 ] 2 

and (2) for the hard stainless steel, 

E [l°gio4>] = - 3. 52 logjo P - 12.20 

°logio4> = [ 2 - 50 + 148 logjo p + 0.033 (log 10 p) 2 J 2 


(150) 

(151) 

(152) 

(153) 


When these results, Eqs. 150 through 153, are substituted into the 
Eq. 29 formula, one gets the contours of equi-probable values of 
puncture-flux cj> versus shell thickness p which are illustrated graphi- 
cally in Figs. 7 and 8 for the hard aluminum alloy and hard stainless 
steel respectively. 


C. THICKNESS OF A FREE WALL VERSUS THE PRODUCT OF 
EXPOSED HEMISPHERICAL AREA AND DURATION FOR 
GIVEN PROBABILITIES OF NO PUNCTURE OF A VEHICLE 
IN A NEAR-EARTH ORBIT 


By substituting the Eq. 138 expression for logj 0 4> in the 
common logarithm of the natural logarithm of Eq. 33 and solving the 
resulting equation explicitly for logj 0 p, one gets 
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logio p = pr 1 [(l - y 5 )/3 + y 5 y 8 + P 3 y 6 ][-log 10 At + log 1 o(-log e R) -yj + 
+ log 10 4.01 + y 1Q + y^ + iy 6 [logio 8. 44 + log 10 (p t /H t )] - 
- y 8 logio Pt + y 6 !ogiocosx 2 +y 3 y 6 + (y 8 - 1/3) [y 4 + 

+ P 5 (P 4 - y s )] • < 154) 

Because log 10 p involves no ratios of random variables, it 
follows by Eqs . 3 and 4 that logio P i n Eq. 154 should stand an even 
better chance of being approximately normally distributed than does 
log 10 4> in Eq. 138, discussed already in Section IV. B. ; especially there 
is more reason for assuming that the mean of the function is the func- 
tion of the means. It will therefore be assumed that Eq. 29 is a 
sufficiently approximate formula for obtaining the antilogarithm of 
Eq. 154 expression for log 10 p. Then, with the same numerical values 
as were used to get Eqs. 139 through 149 from Eq. 138, the mean and 
standard deviation of the common logarithm of the free wall thickness 
p centimers necessary to increase to R the probability of no meteoroid 
puncture of an exposed hemispherical area A square meters during an 
interval of t seconds in orbit near the earth are respectively: 


E [logjo p] = 0. 284 [log 10 At - log 10 (-log e R)] - 0. 288 log 10 p t - 


<r 


logio P 


- 0. 25 log 10 H t - 2. 37 


/ 9 logic P \ 2 , / 9 logio P 

v 8y I "yj \ a* 2 



+ 


v ( B logio p 

L \ ay^ 

i = 3 



2 


+ 


9 logio P 

9y 10 ^10 


2 


/ 9 logio ,P gr f = 0.0974 

\ 9y x Yi) 

( - lo SioP p- x ) Z = 0.0289 
\ 9X 2 2 ' 

/9i£gio_p \ 2 _ o 0072 
\ 9y 3 y 3 J 


(155) 


(156) 


(157) 

(158) 

(159) 
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fl^g UpP r V 

l 9y 4 y* I 

(d logioP V 
\ ay 7 VJ 

fi i S S i PP a V 
\ 9y 10 y W 

/jiogiflje \ 2 
\ 9y g y 5 / 

/9 logxop \ 2 
l 9y fe °V 6 / 

/_ lo gloP q. ^ 2 
\ 9y ^ei 


0.0002 

(160) 

0.0064 

(161) 

0.0100 

(162) 


[0.0129 log 10 At - 0. O1291og 10 (-log e R) - 0. 145] 2 

(163) 

[0. 00382 log 10 At - 0. 00382 log 10 (-log e R) + 

+ 0. 0556 log 10 (Pt/H t ) + 0. 198] 2 (164) 

[0. 0057 logic At - 0. 0057 logic (-log e R) + 

+ 0.086 logic p t + 0.0196] 2 (165) 


These results, Eqs. 154 through 165, applied to each of the 
two metals described by Eqs. 136 and 137, lead to the following 
Eqs. 166 through 169 for formulas for the means and standard devia- 
tions of the common logarithms of the necessary free wall thickness 
p centimeters: (1) for the hard aluminum alloy, 


E [logic p] = 0. 284 [logic At - logic (-log e R)] - 3.03 (166) 

LogioP = {°- 185 - 9* 00229 [logic At - lo gl0 (-log e R)] + 

+ 0. 000213 [logic At - log 10 (-log e R)] z j -2 ( 167 ) 

and ( 2 ) for the hard stainless steel, 

E [l°giop] = 0.284 [logic At - logic (-log e R)] - 3.24 (168) 

°logiop = {0-192 - 0. 00183[log 10 At - logic (-log e R)] + 

+ 0. 00021 3 [logic At - logic (-log e R)] 2 j-f (169) 

When these results, Eqs. 166 through 169, are substituted 
into the Eq. 29 formula, one gets the contours of equi-probable values 
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of necessary equivalent free wall thickness p versus the effective 
area-time exposure product At as illustrated graphically in Figs . 9 
through 11 for the hard aluminum alloy and values of 0. 85, 0.90, and 
0.95 respectively for the puncture-free probability R, and in Figs. 12 
through 14 for the hard stainless steel and values of 0.85, 0.90, and 
0. 95 respectively for R. 

D. VARIATION OF OPERATIONAL, PARAMETERS 

As an illustrative example, consider that one has a vehicle 
with a total surface area of 390 square meters (4200 square feet) of 
0. 318 cm (eighth-inch) free wall of hard aluminum alloy, and one 
wants to know how long the vehicle can remain in a near-earth orbit 
with not less than an even chance that the no-puncture probability 
is as high as 0. 90. 

In Fig. 10, the horizontal line at 0.318 cm, corresponding 
to the eighth-inch wall thickness, intersects the 50% confidence contour 
at the indicated ordinate, corresponding to the 10 7 - 90 square-meter- 
seconds value for At. As the numerical values for meteoroid flux, 
which were established in Section III. A. 4. and used in the derivation 
of the relations which are illustrated in Fig. 10, are based on the 
assumption that only half of the surface area of the vehicle which is 
facing away from the earth will be considered in the flux- area -time 
product, the appropriate value for A in this example is 195 square 
meters, corresponding to the 2100 square feet hemispherical area. 

The permissible exposure time is therefore (10 7 ' 90 / 1 95) seconds, 
which is 4. 8 days. 

E. VARIATION OF DESIGN PARAMETERS 

As a further illustrative example, consider that one has a 
vehicle with a total surface area of 390 square meters (4200 square 
feet), and one wants to know what equivalent free wall thickness of 
hard aluminum alloy will give not less than an even chance that the 
probability of no puncture is not less than 0.90 as the vehicle remains 
in a near-earth orbit for not less than 30 days. The (390) (30) square 
meter days corresponds to the 10 8 ' 70 square-meter -seconds value for 
At, for which an ordinate has been marked on Fig. 10. This ordinate 
crosses the 50% confidence contour at the 0.54cm value for p, 
corresponding to the 0. 21 inch neccessary wall thickness. 
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F . RELATIVE CONTRIBUTIONS TO THE DEGRADATION OF 

CONFIDENCE 

One coincidental but convenient relation between wall thickness 
and confidence can be seen to be common to all of the Figs. 9 through 
14: with both materials for all three values of the no-penetration 
probability, and for all values of the exposure product within the 
interval 10' 1 < At £ 10 u square meter seconds, confidence is increased 
from the 25% value to the 50% value, or from the 50% value to the 75% 

value, accurately by doubling the wall thickness. 

The relative contributions to the uncertainty can best be 
taken as the relative components of the variance of the logarithm of 
the necessary wall thickness, i.e. , by the ratios of each of Eqs. 157 
through 165 to the square of Eq. 156. The numerical values, which 
are indicated in Table I, have been calculated for the specific example 
described in Section IV. E. 


TABLE I 


Parameter 

Source of the Uncertainty 

Magnitude of the 
Uncertainty 


Flux versus Mass 

53.2% 

X 2 

Angular Deviation from the Normal 

15.8 

y 3 

Closing Velocity 

3. 9 

y 4 

Mas s -Dependenc e of Density 

0. 1 

y 5 

Mass-Dependence of Density 

0. 2 

y 6 

Penetration Exponent of Velocity 

10. 8 

y 7 

Penetration Coefficient 

3.5 

y 8 

Penetration versus Target Density 

7. 0 

y 

10 

Target Thickness Factor 

5.5 
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The last four lines in Table I indicate the relative components 
(of the variance of the logarithm of the necessary wall thickness) which 
are due to our inability to predict more accurately what will happen in 
a specified hypervelocity impact; and the sum total of these four 
components is only 26.8%. But if these four components could be 
entirely eliminated, then the probable error in the logarithm of the 
necessary wall thickness would only be reduced to 85. 6% of its present 
value; i. e. , ( 1. 000 - 0.268)^ = 0.856. 

G. AMELIORATING CONSIDERATIONS 

Lest one might be apprehensive of the considerable free wall 
thicknesses which are inferred in Figs. 5 through 14, it is appropriate 
to note that Nysmith and Summers [33] report: . .two-sheet structures 

are probably the most efficient multiple -sheet structures. The most 
substantial gain in penetration resistance is achieved by filling the void 
between the sheets with a glass-wool filler. For the one case, investi- 
gated, namely, a two-sheet structure with a sheet spacing of one inch, 
the penetration resistance of the structure with the glass-wool filler is 
about twice as great as that of the structure without the filler material 
and about 4.4 times greater than that of a single sheet of material of 
the same total-sheet thickness. M 


SECTION V. CONCLUSIONS 


It can be seen from the numerical results, which are graphically 
illustrated in Figs. 5 through 14, that if a vehicle which is very large, 
and which also has thin walls, remains very long in orbit near the 
earth, it may be punctured by a meteoroid. 

For a vehicle of given size, given time in orbit, and given 
no-puncture probability, it is necessary to increase the effective free- 
wall thickness by a factor of 2. 00 to change the confidence from 50% 
to 7 5%. If hypervelocity puncture effects on given materials were 
completely predictable for projectiles of specified density and velocity, 
then the above factor would be 1. 81 instead of 2. 00. This would help 
appreciably, but not so much as would a reduction in the uncertainty of 
the mean-time -between-hits per unit of surface area for meteoroids of 
mass equal to or greater than stated values. 

The reader should be aware that there still exists some 
controversy over the available meteoroid flux, density, velocity 
and angular distribution statistics. Therefore, for any design 
studies the latest accepted values should be employed. 
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l°g 10 m 

FIG. 1. MASS DEPENDENCE OF THE FLUX OF METEOROIDS ON A SPHERICAL SURFACE 
IN A NEAR EARTH ORBIT. MASS >m (gm). AVERAGE NUMBER PER SECOND 
PER SQUARE METER AT 25, 50, AND 75% CONFIDENCE. 











cm arp at Confidence 



m (gm) 

FIG. 5. MASS DEPENDENCE OF THE THICKNESS OF AN EMPTY SHELL 
OF HARD ALUMINUM ALLOY JUST PUNCTURABLE BY A 
METEOROID OF MASS m 
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& p at Confidence 



m (gm) 


FIG. 6. MASS DEPENDENCE OF THE THICKNESS OF AN EMPTY SHELL 
OF HARD STAINLESS STEEL JUST PUNCTURABLE BY A 
METEOROID OF MASS m 
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centimeters 



10" 12 10 -8 10 -4 10 0 

sec" 1 meter" 2 * <j> at Confidence .25, .50, and .75 


FIG. 7. METEOROID PUNCTURE FLUX FOR A FREE WALL OF HARD 

ALUM I N UM ALLOY FOR A VEH I CLE I N A NEAR-EARTH OR B I T 
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centimeters 



FIG 8. METEOROID PUNCTURE FLUX FOR A FREE WALL OF HARD 

STAINLESS STEEL FOR A VEHICLE IN A NEAR-EARTH ORBIT 
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cm * p at Confidence . 25, . 50, and . 75 



FIG. 9. NECESSARY THICKNESS OF HARD ALUMINUM ALLOY FOR 
.85 PROBABILITY OF NO PUNCTURE OF A SPHERICAL 
VEHICLE IN A NEAR-EARTH ORBIT VERSUS EXPOSURE At 


6k 




a p at Confidence .25, .50, and .75 











cm kpat Confidence . 25, . 50, and , 75 



At, Square Meter Seconds of Exposed Hemisphere 


FIG. 11. NECESSARY THICKNESS OF HARD ALUMINUM ALLOY FOR 
.95 PROBABILITY OF NO PUNCTURE OF A SPHERICAL 
VEHICLE IN A NEAR EARTH ORBIT VERSUS EXPOSURE At 
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at Confidence .25, .50, and .75 



10' 1 10 1 10 3 10 5 10 7 10 9 10 11 
At, Square Meter Seconds of Exposed Hemishpere 


FIG. 12. NECESSARY THICKNESS OF HARD STAINlfSS STEEL FOR 
.85 PROBABILITY OF NO PUNCTURE OF A SPHERICAL 
VEHICLE IN A NEAR-EARTH ORBIT VERSUS EXPOSURE At 



p at Confidence . 25, . 50, and . 75 










at Confidence . 25, . 50, and , 75 













ANTIHELION 

A-Al'80 1 


A-A s - 270° 6i 


Time of I Transit 
T 


EARTH 



ANTAPEX 


c 8h '6h ij A- A s =90° 


A-A|0' 


FIG. 15. POLAR DIAGRAM DRAWN IN THE PLANE OF THE 
EARTH'S ORBIT WHICH SHOWS THE APPARENT 
NUMBER OF METEOR RADIANTS DETECTED PER 
UNIT ANGLE, PER SECOND (REPRODUCTION OF 
FIGURE 7 OF REFERENCE 51) 





ANTIHEUON 



FIG 16. POLAR DIAGRAM DRAWN IN THE PLANE OF THE 
EARTH'S ORBIT WHICH SHOWS THE NUMBER OF 
METEORS PER UNIT ANGLE WHICH CROSS THE 
EARTH'S ORBIT PER SECOND (REPRODUCTION 
OF FIGURE 8 OF REFERENCE 51) 
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